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ABSTRACT 


It is shown by direct calculation to fourth order in the coupling constants 
that all divergences, except those arising from the scattering of mesons by 
mesons, can be removed from the combined interaction of charged spinless 
mesons, nucleons (scalar interaction) and the electromagnetic field by mass 
and * charge ’ renormalization. The remaining divergence can be cancelled 
to this order by assuming a direct interaction between mesons with a 
suitably chosen coupling constant. 


§ 1. INTRODUCTION. 


Dyson (1949) has put Feynman’s (1949) graphical technique for calculating 
the S-matrix elements in electrodynamics into the Schwinger (1948) 
Tomonaga (1946) form and has shown that the consistent application of 
mass and charge renormalization yields a finite theory to any order in the 
coupling constant. The analogous treatment of the scalar interaction of 
spinless mesons with nucleons has been given by the author (Matthews 
1950 a) and also a summary of the graphical method as applied to charged 
spinless mesons in the electromagnetic field (1950 b). This latter inter- 
action has been treated independently and in much greater detail by 
Rohrlich (1950). Here the extent to which renormalization is effective in 
removing the divergences from the combined interaction of neutrons, 
protons, charged spinless mesons (scalar interaction) and the electromagnet 
field is considered by direct calculation to fourth order in the coupling 
constants e and f. It is shown that for the removal of all divergencies, 
except those from the scattering of mesons by mesons, by mass and 
“ charge’ renormalization certain conditions between the infinite constants, 
which determine the renormalized charge, have to be satisfied. It is 
checked that these conditions are satisfied to second order. The remaining 
divergencies from the scattering of mesons by mesons can certainly be 
removed, to this order, by the introduction of a direct interaction between 
mesons with a suitably chosen constant thus leading to a finite S-matrix. 
(A brief report of these results has been published elsewhere (Matthews 
1950 ¢).) 

One of the conditions referred to above occurs in the interaction of spin- 
less mesons with the electromagnetic field. Since the completion of these 
calculations, more general proofs of it have been given by Dyson (1950) and 
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by Rohrlich (1950), based on gauge invariance. A similarly general 
solution has since been given by Salam (1950) of the present problem by a 
development of the work of Ward (1950). The problem in its complete 
generality is much complicated by the overlapping of divergent parts 
similar in principle to the ‘“‘a’’ and “‘b” divergencies in self energy parts 
discussed by Dyson (1949, § VII) and no proof has yet been given, though 
there seems little doubt that this scheme does in fact lead to a finite 
S-matrix up to any order. 


§2 PRIMITIVE DIVERGENTS. 


The Hamiltonian in the interaction representation of the combined 
interaction for pseudo-scalar mesons is 


H=2H,,.,.(==)5/2 5.3) q;uek dee ween 
where 
db  ad* ie \2 
H,= (=) A(t 52-54) —(F) MAb bon 2) 
Hy=iflys(7_$* +7. dp + S692 ernest «cE 
H,=tepy, TppA,,. sk ows: ocd ces scone ieee ea etal eo 


7, operating on a neutron gives a proton and on a proton gives zero. 
7_ is the same operator with the neutron and proton reversed. Tp gives a 
proton when operating on a proton and zero when operating on a neutron. 
The terms required for mass renormalization have not been included. A 
term in H, explicitly dependent on the surface direction has also been 
omitted. This can be ignored without error if the corresponding term in 
the P-bracket of two derivatives of ¢ is also ignored (Matthews 1949.) 
Any scattering effect in which both the real and virtual processes are specified 
can be represented by a graph. Dotted lines denote photons, wavy lines with 
double arrow denote mesons, full lines with double arrows denote protons 
and full lines with single arrows denote neutrons. Thus double arrows 
signify the movement of charge and for charge to be conserved the double 
arrowed lines (either full or wavy) must form continuous loops (open or 
closed). Single arrows determine the motion of neutrons (as opposed to 
anti-neutrons) and for consistency the nucleon lines must form continuous 
loops of either single or double arrows. At any vertex with a photon line 
entering the other two lines must be double arrowed. At a meson- 
nucleon vertex the nucleon lines must be one double and one single arrow. 
If these rules are observed the 7-factors can be neglected. The rules for 
obtaining the corresponding S-matrix element can be summarized as 
follows (see Dyson (1949), Rohrlich (1950) and Matthews (1950 a and b)) : 


(i.) Each internal photon line give a factor 
hed,,y(27)-8 [ Dp(p)d*p. 

(ii.) Each integral meson line gives a factor 
he(2r)-* | Ay (p)d4p, 


‘Be 
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(iii.) Each internal nucleon line gives a factor 
(27r)~8 [Sp(p)d4p. 
(iv.) Each meson—photon 3-vertex gives a factor 


te(he*(27)*(p,+-p,)3(P—p' +9), 
where p and p’ are the four-vectors of the meson lines entering and leaving 
the vertex, respectively. (The factor (p,.+p,) arises from the derivatives 
in H,.) 
(v.) Each meson—photon 4-vertex gives a factor 


—te*(he)?(277)48,,,8(p—p'+9—q’). 
(vi.) Each proton—photon vertex gives a factor 
e(fic)~1(27)4y,,. 
(vii.) Each meson-nucleon vertex gives a factor 
(fie) "(2m)4y5. 
_ (viii.) Each external photon line gives A,(k), each external meson line 
entering or leaving the graph gives ¢(k) or ¢*(k), respectively. Each 
external nucleon line entering or leaving the graph gives %(k) or i(k), 
respectively. 


(ix.) There is a further numerical factor of (—1)"~"-1, where n is the 
number of nucleon vertices and / is the number of closed nucleon loops. 
Finally, the whole integral is multiplied by the number of ways the 
operators can be paired off by interchanging the roles of the two photon 
operators at the 4-vertices. 

The condition for primitive divergent graphs can be derived by Dyson’s 
method (1949) and is 

(3/2) ele ls eee soe = (3) 
where E,, E,,, and E,, are the numbers of external nucleon, meson and 
photon lines, respectively. Since the nucleon factors appear in H in pairs, 
E,, in (5) must be either two or zero. If E,, is two then for (5) to be satisfied, 
either E,, or E,, must be zero. Since the charge bearing lines in any graph 
form continuous loops, the number of external charge bearing (double arrowed) 
lines in any part must be event. These considerations together show that 
the primitive divergents for the three field mixture, in terms of their external 
lines, are the same as for the separate pairs of interacting fields. They 
may, however, have different internal structure. Some of the new graphs 
are illustrated in fig. 1. (It should be noted that these all satisfy the 
charge conservation and consistency conditions given above). 


ee eee ee ee eee 

+In the author’s treatment of the meson—nucleon interaction? (1950 a), 
it was wrongly stated that for neutral mesons the divergences from parts 
with three external meson lines cancel by a type of Furry theorem (Furry 
1937). Genuine divergences may arise from this type of graph for neutral 
mesons, but are definitely excluded for charged mesons by the charge 
conservation rule, 
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§3. RENORMALIZATION CONDITIONS. 


Considerations are now restricted to fourth order in the coupling 
constants e and f. To calculate the probability of any process one draws 
all the irreducible graphs of the process to this order. It should be 
remembered that C-parts (with two external meson and two external 
photon lines) are divergent. Any such part (except that with a single 
internal meson line; fig. 2 (a)) is replacable by a 4-vertex and counted as 
reducible. To obtain the corrections, after mass renormalization, due to 
the insertions of self-energy, vertex and C-parts, D, and 4, must be re- 
placed by Dj, and 4; defined as 


Dy(p)= (145A) Da()+DR(P), ss 
Ayly)= (1+ C+ EC) Aylp}+ ARC), 
Fig. 1 
(a) (b) 
(c) 
(d) eS 


where o=e?/4rhic and B=f?/4nfic. A, © and C’ are infinite constants 
defined analogously to those occurring in Dyson’s treatment of electro- 
dynamics. D) and 4p) are finite and of second order in the 
coupling constants. The 8,() factors for proton and neutron lines must 

be replaced, respectively, by 


SHo)= (14+ EB+ EB) S(H)+SKP(r), (8) 


Si(p)= (1+ 5B) Sp) +892 (p). 2,714 201 ere O) 


of the Meson—Photon—Nucleon Interaction 225 


.Also the vertex factors ¥5> Yw (Pup) and 6, must be replaced 


y aN a Me B Wi yy i? 
LS(p, p’) (1 a, PA FP’) vo bAB (m,p Vector Congres aera LO) 
Fa g\ = B , y / 
LAD; Pp )= (1+ opie sh) nbae (p, Pp ), ORO Oe Cie 0 (11) 
®,(p, p= (1+ = N+ LN’) (p,+p)4+ 0% (p, p' 
WNL? o- ay on (Put Pu) + Or (Dp, p ), B Te (12) 
GM ges (pak Pe, 
Be?) a RRO) 8 4 OO) (p, pg), 2 9. (13) 


P, P’, L, L’, ete. are infinite constants and the final terms in these equations 
are all finite. The constants B and L have been calculated explicitly by 
Karplus and Kroll (1950). The ‘infinite parts” of the constants C, N and 
R have been calculated independently by Rohrlich (1950). Also external 
line factors must be replaced by (I +«/47rA)A(k), for a photon, 
(l--a/4rC-+-B/47C’)4(k) for a meson, (1+-«/47B-+8/47B’ yb for a proton and 
(1-++8/47B’) for a neutron. 


Fig. 2. 
(a) (b) 


The idea of charge renormalization is to express any infinite term occur- 
ring in the calculation of a matrix element as a constant infinite multiple 
of a finite term of lower order and then to show that these infinite multiples 
can be absorbed by defining the observed coupling constants e, and f, as 
double power series ine and f. To see whether this can be done consistently 
in this case, consider the scattering of photons by mesons. The only 
irreducible graphs are shown in fig. 2. For renormalization to be 
effective, the infinite multiples of these graphs arising from the substitutions 
(7) and (12) into fig. 2 (a) and (13) into fig. 2(b) must be the same, since 
each graph is of order e?. There are also infinite multiples arising from 
substitutions in external lines, but since these are the same for both graphs 
they need not be considered. The required conditions are 


CNC ae a. aml) 
ON" Cee Rie joc o. | Suen 


Condition (14) allows for insertions involving mesons and photons only and 
has been checked by Rohrlich (1950) and by the author (1950 b). Condition 
(15) arises from the additional parts due to virtual nucleons. (The term 
2N’ comes from the insertion of fig. 1 (a) in the two vertices of fig 2 (a), 
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C’ from a megon self energy part due to virtual nucleons in the internal line 
of fig. 2 (a), and R’ from the insertion of fig. 1 (b) in the 4-vertex of fig. 2 (b). 
This condition has also been checked. ‘The required charge renormali- 
zation is 


ane[1+¢(n+c+5) +0409], ye ES 


since each factor e is associated with a 3-vertex a meson line and half 
a photon line. As has just been verified, this same renormalization absorbs 
the factors R and R’ from the C-parts. 

If processes involving real protons and photons are now considered, 
other infinities will arise from the substitutions (6), (8) and (11). These 
are removed by the charge renormalization 


ane 1+ ¢(BHL+ 5) +5 B41) |. Oe ae 


For consistency this must be identical with (16). That is, it is required 
that 
BL LENS OM ea es 
B/+L’/=N’+C’. Sy nth Ai aae el © eee 
In fact it is found that 
B=—L, Bo=—Li, N=—C=R, 7 N=—CH]R. eg) 
This means that the whole effective charge renormalization comes from 
photon self energy parts which are the same for both meson—photon and 
proton—photon effects, thus leading to the same charge renormalization 
for both. 

For meson-nucleon effects the infinities arising from the substitutions 

(7), (8), (9) and (10) can be removed by the renormalization 
, 
f=t|1+3 - (5 =F = 5 +?) + 5 (B45 +P) Seine caer 6-4, 
This can be made independently of (16). 

A further possible source of div ee is from graphs with one photon 
and two neutron external lines (fig. 1(d) and (e)). Since no (finite) 
single vertex of this type exists no renormalization is possible. 

However, though each graph diverges separately, the total is finite and 
no difficulty arises. These graphs lead to the finite neutron magnetic 
moment which has been discussed in this framework by Case (1949)). 

The only remaining divergencies are those from the scattering of mesons 
by mesons (M-parts), with four external meson lines. This can occur 
either through virtual nucleons, (f4), or virtual photons, (e4), and in,each 
case gives a genuine divergence. To this order these divergences can be 
cancelled by a suitable choice of the constant in the direct interaction 
dAd**4?. This device has been considered more generally by the author 
(1950a) for the meson-nucleon interaction and by Rohrlich (1950) for the 
meson—photon interaction. The former treatment was very incomplete 
because secondary infinities from parts such as fig. 3 were not considered. 
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Rohrlich’s development is more thorough but there are difficulties when 
such parts appear inside self energy parts which are not dealt with in 
' detail. This point will be discussed further in a forthcoming paper by 
Mr. A. Salam: 

The whole of the discussion is algo applicable to scalar mesons. (It is 
only necessary to replace (Ys)ag by —18,,-) 


Fig. 


3. 
“ae eG 
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The calculation of the infinite constants was performed by the methods 
due to Feynman (1949) as employed by Karplus and Kroll (1950). The 
“infinite parts”’ of all the constants can be expressed in terms of 

eel d*k 
Dae p +2 
and. are 
A=—3l P=i1,.P’=0, 
N=—C=R=—I, B=—L=4I, 
Ni=—C R=] Bo=—L'=4HT. 
(It is readily seen that the single and double arrows cannot be inserted 
consistently in the graph of this order corresponding to P’. This is not 
true for higher orders.) 
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to Non-linear Differential Equations. 
By A. C. Sim, 


Standard Telecommunication Laboratories Ltd.* 


[Received January 11, 1951.] 


SUMMARY. 


The formule for algebraic reversion are extended to revert a class of 
non-integral power series, and also generalized to revert series in which 
the coefficients contain operators. These formulz are shown to have a 
simple and powerful application to non-linear differential equations, which 
they expand into an infinite sequence of linear equations. 


§1. IvrrRopuction. 


THE studies of semiconductors and saturable iron-cored coils are two of 
many important fields in modern physics and engineering which lead to 
non-linear problems of practical importance. Such problems occur also, 
for example, in aerodynamic and navigational theories, in the analysis 
of electronic circuits, and in nuclear physics. A comprehensive survey 
of the theory of non-linear mechanics has been published by Minorsky 
(1947), from which it is clear that the number of methods which are not 
restricted to one type of problem is very small, and most of the remainder 
are concerned with non-linear oscillatory equations. The method of 
reversion described herein is applicable to a wide variety of equations, 
of which some appear otherwise intractable, provided the equations are 
expansible in powers of the unknown quantity. 

The existing reversion formule are applicable only to algebraic 
equations, and are used extensively for the calculation of roots. The 
problem may be defined as follows. Given the series 


A= 2G, .Y'; « rohan te, fetes es a 
the coefficients A, in the reverted series expansion 
Ya 2 AN een an eee) 
r=1 


have to be found. Van Orstrand (1910) has published a very full account 
of such formule and their applications. In addition, a special type of 
series occurring in the theory of Bessel functions and having the form 

{Yat YAO ete eee apne) 
has been reverted by Bickley and Miller (1943). 


* Communicated by Professor D, Gabor. 
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If the series (1) and (2) converge * it is clear that (2) is a solution of 
the equation (1), but not necessarily the only solution. In fact it is the 
solution which vanishes as approaches zero, and is, in general, the smallest 
real solution. If equation (1) possesses no real solution (2) will diverge. 

The reversion formule may be obtained in two ways, of which the most 
obvious is the substitution of (2) into (1). The quantities A, are then 
found by equating the coefficients of \” on each side of the equation. A 
more elegant approach, due to Harkness and Morley (1898), involves 
differentiating (2) with respect to y and multiplying each side byeas 
This gives 


A= dAfdy{A,AT+....+7rA A+. 2.0}. 6. Oe (4) 
It is then argued that the relation d” .dA/dy=(n+1)-. dAv+4/dy is true 
for all n except n=—1; and since A”+! is a power series in y, its 


derivative cannot contain a term proportional to A~!. dA/dy, because this 
can only arise from the differentiation of log A. It follows therefore, 
from (4), that the coefficient of y~! in the power series expansion of A~” 
is rA,, and hence 


rA,—coetticient of y1in {a,y+agy?+....} 7 . . . (5) 


It is clear that the evaluation of (5) will assume that the coefficients 
a, commute, and this method therefore fails in general when operators 
are involved. _ 


§2. REVERSION. 
The series (1) and (3) described above are special cases of the non- 
integral power series 


aggre ay isn tie. a fe SUPT aa hes }=A¢(a). ; (6) 


This series is the most general one that can be reverted. The coefficients 
a, may be any functions of x, d/dx, or any other operator, but they must 
be independent of y. The real quantity « is quite arbitrary; is an 
integer—positive, negative, or zero , and A is a constant. 

In the reversion of (6) it is most convenient to consider y+" as the 
unknown quantity. Inspection shows that the reverted series will be in 
powers of y” where m=«/(1-+nx), so that the solution may be assumed 


to have the form 
ye) Ned ha Oe eee Se i ae (7) 
r=0 


In determining A, the extended multinomial theorem is required. 
This may be written in the form 


rs eee LLY) ents nucle 
fuk+uke+...+u,ket... j= 2 2 TiL+4) Pa! Fond. 


r=0 p 


usher, 


(8) 


* Harkness and Morley (1893) have shown that the reverted series (2) 
converges in the same domain as the original series (1). 
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where the exponents g and p, assume all values which satisfy the two 
relations 


(Potpst...-+p,)=("—9); 
(pet 2pst..- + (r—Lp,)=r. 


It should be noted that p, and (v—q) are positive integers, and when v 
is negative, the ratio [(1+v)/I'(1+q) in (8) should be replaced by 


(—1)"-9. P(—g/P(—») 
By direct substitution of (7) into (6) using (8), the reversion coefficients 
are found to be given by 
ty -Ag=$(2), 
dp. Aga—a,, Ags 
dy. Ay=— {a,.(1+m)A”A,+a,. At?” }, 
Ay. Ag= — {a, .(1+m)[mAM-YAZ/2 ! +-AGAd] 
+. (1+2m) AA, +a,.A0t3™}, 
dy -Ag=— {a, . (1-+m)[m(m—1)AM—A3/3 ! +-mAS"— YA, A,+Ap'As ] 
+a. (1+2m)[2mA2™—DAZ/2 | + Ad™A,] 
+a,.(1+3m)Ai"A,+a,. A974}, 


and the general term 


ee I(2+ms) yy», ASS 
Goh ge a Corea cama aes Pat? te Dy 
where (PitPoet..... +Dp,—s) =(1-+-ms— pp) 
and (Pp, +2po+....+(r—s)p,_.)=(r—S). 


Note that (1--ms—p ) is always a positive integer. 
In the general term of (9) when (1+-ms) is negative the ratio 


I(2+ms)/I'(1+ 9) 

should be replaced by (—1)4+™s+»)). '(—p,)/I'(—1—ms). 
In general, the formule given by (9) represent a sequence of linear 
ordinary differential, or integro-differential equations determining the 
functions A,. The equation ay). A >—4(x) is referred to as the generating 
equation and its solution A, is termed the generating function since the 
complete solution is expressed in terms of it. The generating function 
is a first approximation to the true solution and the generating equation 
is a linear approximation of the original non-linear equation. 

It may be noted that if the coefficients a, commute a simplified formula 
can be obtained by the method of Harkness and Morley. It is then 
convenient to change the notation, writing the series in the form 


(yt Gy2-bo ew cpa,et ie es JREASS Wee tt we oe (10) 


ae —eregrnan 
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in which B has replaced (n-+-a-1) and z is written for y*. The function 
#(x) is now included in A. Any arbitrary power of z may be found in this 
case so that a solution may be assumed in the form 


=(Naym SBA | ad) 
r=0 
es ee 
"sin 058: Rp a im 
Bo=1, By=—sm.ay ta, 
B,=sm{(l+2+s. m)ag O+Mat/2! —az it2mMq,), 
B,=—sm{(14+3-+8.m)(2+3-+8.m)ay24+ma3/3 | 
—(14+3-+8.m)a> BRS Mp et bead aa 
B,=sm{(1+4+s.m)(2+4+s8.m)(3+4-+s. mag “AtMatl4 | 
—(14+4+8.m)(2+4+5.m)az @+*™ea2a,/2! 
+(1+4+8.m)aj2+2™(a,a5-+02/2 N—ag Ut *™a,}, 
and the general term 
pee ee (Pot em) ac™, Gite. OFF 
p T(s+r.m) Pyles. Ppl 
where (Pi tpeot....+p,), =(Po—7TmM, 
(p,+2p.+....+rp,) =r. 
In this case when (s+-r)m becomes negative the quantity 


(—1)(2-7”). P(pp-+-sm)/I'(s-+r .m) 


must be replaced by the ratio ['(1—s-+r.m)/I’(1—pp—sm). 

Whilst Van Orstrand has dealt adequately with the application of 
reversion formule to equations not involving operators, it may be noted 
that the extended form of equation given by (10) offers alternative forms 
of solution, and will also yield results when applied to some equations 
containing fractional powers. For example, the equation 


(x3 +-c-+-0,2) = Cy 


reverts to a series in powers of cy. It may be, however, that cy is very 
large and c, is the smallest coefficient. Division of the equation by x 
gives (a?+c,4—cox-1)=—c¢, which may also be reverted, producing a 
more convergent series. It is also possible, using (12) to revert the 
equation (c,x"+c,x°)=cy for any values of @ and b provided they are 


real. 
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§3. Non-Linear EQuatTIons. 


To apply the general reversion formule (9) to a non-linear differential 
equation, the equation must first be expanded in an ascending or descend- 
ing power series of the unknown function y—the coefficients being expressed 
in terms of « and D=d/dz, and arranged so that they contain no negative 
powers of D. Ifa general solution is required, the generating coefficient 
a, Which is the first one in the series, must be of the highest degree in 
D so that the generating equation has the same order as the equation to 
be solved. 

The set of equations given by the reversion formulz (9) are all similar 
in form to the generating equation, although the functions of a on the 
right-hand side become progressively more complicated as 7 increases, 
Each equation, therefore, will give rise to the same number of arbitrary 
constants. For a complete solution of the equation to be reverted, a 
general solution of the generating equation must be found, whilst only 
particular solutions of the remaining equations for A, are required. The 
choice of these particular solutions affects not only the form of the series 
obtained, but also the ease with which the arbitrary constants are deter- 
mined from the boundary conditions of the problem, and the convergence 
and complexity of the reverted series. 

If the reverted solution of an equation does not contain sufficient 
arbitrary constants, or does not converge in a convenient range, the form 
of solution may frequently be improved by a change of variable. In all 
applications, therefore, every conceivable form of the expanded equation 
should be considered, to decide which one is most advantageous. 
Convergence of the reverted solution must be examined in each individual 
case, since a general discussion does not appear to be possible. Because 
of the complexity in deriving a large number of terms of the reverted 
series, it is usually easiest to check its accuracy and convergence by 
substituting the few terms obtained into the equation, and examining the 
remainder which results. Even this procedure can be very formidable, 
in which case a numerical substitution into the equation affords the only 
practical verification. 


§4. EQuaATIONS witH INITIAL BOUNDARY CONDITIONS. 

It is generally the case that the arbitrary constants in the solution of a 
non-linear equation are so woven within it that their determination becomes 
a difficult task. In this section it will be shown that when the generating 
equation is linear with constant coefficients, as is frequently the case in 
practise, a form of solution can always be found by reversion which makes 
the calculation of arbitrary constants from initial conditions as easy as in 
the case of a linear equation. This is achieved by choosing the particular 
solutions for the functions A, so that they all vanish (and also their first 
n—1 derivatives where n is the order of the equation), when the inde- 
pendent variable is zero. The constants can then be rigorously calculated 
by regarding the generating function as the complete solution for the 
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purpose. The particular form of the generating function thus obtained 
can then be used for the derivation of the subsequent terms in the series 
until sufficient accuracy is obtained from them. 

It is now assumed that the generating coefficient @ 18 a polynomial in 
powers of D=d/dx with constant coefficients. If the degree of this 
polynomial (equal to the order of the non-linear equation), is », the 
inverse of a may be expanded into partial fractions. Let these partial 
fractions have the form 


, = Vo t rs 8 > 
ee 2 anemia eens 
gmt Dt) | snr gui (+B) ee 
in which —a,, —%,...., are the non-repeating zeros of a», and 
—P;, —B,,...., are t roots which are repeated r,, r,,...., ¢, times. 


The coefficients y, and 5, ,, are constants depending on the polynomial 
coefficients which are not defined here. The quantity R is the difference 
between the degree of the polynomial and the number of non-repeating 
roots, 7.e. R=2y7,. 

The reversion functions A, are all defined by the operation of the 
collection of operators given by (13) on a function of 2 which is known 
or can be explicitly determined. Let f(x) represent this function for a 
particular A,. It is desired to find a particular solution to this operation 
which will vanish with its first »—1 derivatives when «x is zero, then 


Da, 70) 1s required tor u—0, 1,...,2—1. 
With the aid of the contracting notation 
v (—lI)¢" .u! 7 
O,(bs—-0,) = 3 ae Del Mie M/S) (14) 


it can be found by means of simple algebraic manceuvres that 


5 peer =8, 5 6B) DH e+6,—B YF feoeh. . (5 
pei 8,)2 ee : at Ete 


Thus in the determination of A, the only operations that need to be 
considered are D2%.f(x) and (D+8,)~¢.f(«). The first of these is simply 
a repeated differentiation, whilst the second is readily interpreted using 
the operator shifting theorem. A solution which vanishes when « is zero 
‘can thus be written 
% nie 
(D+B,)-*. fae" | de| tg le (¢ times) a pete f(x). 00... (LG) 
0 0 . 

Since a, is defined to be a polynomial of degree n in D, if ay ‘is expanded 
in powers of D~! the series obtained will commence with a term in D™”. 
But, by expanding the right-hand side of (13) in powers of D~! some 
terms containing lower powers are obtained. By equating the two 
expansions it follows that the coefficients of D~* in the expansion of (13) 
for k=1, ..., n—1, are zero. In the notation of (14) this result can be 


found to be 
n—R t a ai ; 
{ = shyt 2 (-1)5,08,—B.} =0 k=0,1,...,(m—2). (17) 
\ aoe. 
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It now follows from (13) and (15) that 


n—-1n—R 
Ppa Agee { EE (—1fkaby, De, f(v) 


k=0 q=1 


+ 25,5 @,-P)D-. fla)} 
z=0 


s=1 k=0 


u=0, 1. 5 (1) oe ee 


Finally, from (17) and (18) it is seen that if a solution for A, is obtained 
in this special case that vanishes when a0, then its first (m—1) 
derivatives will also vanish when z=—0, 7. e. 

If A,=0 when «=0, then 


(2*A,/d8") p= 0 for U1, 2, ee ee ee ed 


when dy 1s a polynomial of nth degree in D with constant coefficients. 


§5. A SPECIAL EQUATION. . 


Apart from the general power series (6), it is possible to revert the 
equation : 


P.y={e)tf@)} . 6 2... -(20) 


when f(y) is any function of x and y, but ¢(x) is independent of y as before. 
P is any operator. Note that it is not possible to extract a common 
factor A in this case. This form of reversion may be useful when the 
equation is not expansible in powers of y, since the procedure does not 
demand this. To solve the equation a solution may be assumed to exist 
in the form 


y={P-1. d(z)+F(z)} . 5 2 2... QD) 


and this may be substituted into (20). A power series in F(x) is then 
obtained by expanding f(y) using Taylor’s theorem. The resulting series 
can then be reverted by means of the formulz already derived in § 2 of 
this paper. The final reverted solution is then 


yas 2 OP tee ae a 
r=0 
in which 
P.Cy=¢(zx), 
(P—f) . C,=f(Co), 
(P—f'). Ca=fl > Ci/2 1, 
(P=27)) 7. Case hart Aes Goeey cae C3/3 !} ete. fibh itupalea) 


with the abbreviation f” = E fw) | i 
ar y=Cy 
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A similar set of formule may be derived for the case when (20) is 
without operators, and it is easy to verify that the result is equivalent 
to that obtained by Lagrangian expansion (Whittaker and Watson 1950) 


§6. ILLUSTRATIONS. 

(1) Semiconducting materials are sometimes used for shunting switches 
to reduce transients in the breaking of inductive circuits. In the circuit 
of fig. 1, an inductive coil is represented by resistance R and inductance L 
in series, and it is supplied from a battery of voltage E. The switch S 
is shunted by a non-linear resistor possessing the characteristic V—KI-, 
where V is the voltage across it, I the current through it, and K and c 
are constants. The equation satisfied by the current may be written 


Tet Dieler i egithes eas |. | (24) 


where the constant a represents the ratio K/R, and the operator D is 
equivalent to L/R .d/dt. I) is the initial value of I, equal to E/R. 


Biget: 


This equation has the special form discussed in § 4. The generating 

equation is (D+1) . I=I, so that the generating function is 
I,(1—exp [—R¢/L])-+A exp [—R¢/L]. | 

where A is a constant. Regarding this as the complete solution for the 

purpose of determining A, it is found that the exponential terms vanish, 

so that in this case the generating function is simply ly ee us and 

equation (9), the next term in the reverted series is —(D+1)-1. aIj, and 

the series is found to be 

T= {1,—al%(1— exp[—Re/L])+ ca? Ife (1 — exp[—R#/L](1 +R¢/L)) . Ae 

| ”) 

As in the case of algebraic reversion, if the unknown quanuty Tihs 

replaced by a known approximate solution plus a new dependent pages 

the reverted series will then converge more rapidly. In this problem i 

is known that the initial value of I is Ip, so that the ges of 

(I’--I,) for J should improve convergence. This transforms the equation to 


f(D AST V's a(T ely hee ee ene Oy 
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Defining parameters p and q by 
p=[l1+acl¢-Y] g=ac(e—1)1f-/2 
and expanding the second term in (26) by the binomial theorem puts the 
equation in a suitable form for reversion 
{(D+p)I'+qI?+...}=—alj. . - . . . (27) 
Reversion of this series produces a series more rapidly convergent when 
a and c are large, and no less convergent when they are small 


I= {I,—ap-11g(1 —exp [—pRt/L]) 
—a®p-8q1‘(1—exp [—2pRt/L]—2Rt/L . exp [—pRt/L]) ..-}. (28) 
Fig. 2. 


Differential 
Analyser, 


Equation (29). 


0.05 0.10 0.15 0.20 0.25 0. 30 
Rt/L 
This problem has been analysed by Fairweather and Ingham (1941) 
using a differential analyser, and it is interesting to compare equation (28) 
with their results. As a particular numerical case, let c—0-2, K—240, 
R=50 and E=50. Then I,=1, a=4:8, p=1-96, and q=—0-384. 
Equation (28) then assumes the form 
= {1—2-45(1—exp [—1-96Rt/L]) 
— 1-182(1—3-92 exp [—1-96Rt/L]—exp [—3-92R#/L])...}.. (29) 
In fig. 2 a plot of this expression is compared with the curve obtained on 
the differential analyser. 
(2) The functional equation* f(w+1)=e/ may be written in the form 
e?. f(x)=e/, and with y=f(«) this may be expanded into the series 


{(e?—1) ge y2t—y/al. teed we BO) 


a a er ee ee ee eee 
* This problem was communicated by Professor K, W, Miller of the Armour 
Research Foundation, Chicago, Ill, 
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A solution of this equation is required that vanishes when v—0. Using 
the expansion 


(e°—1)-47=D-1{1—D/2-+ D412... —(—1)"D® , B, /(an) 1}, 


where B,, is the nth Bernoulli number, the generating function is seen 
to be Ap=(e?—1)-1.. 1=(x+a) where a is a constant, and in this case 
a=0. Proceeding with the reversion leads to 


Y=} (x)= {x+a(1—2x)[(1—2x)/12+ (2—13¢@-+1722 8x°)/240...J}. . (31) 


(3) The performance of an anharmonic oscillator satisfies the equation 

Desa) yt by 0, ey, 32) 

For a general solution the generating function is A cos (ax-+-c) where 

A and ¢ are constants. Reversion then produces the series 

y= {A cos (ax-+c)—bA?/2a7[1 —cos (ax)-+cos (ax+2c)/3 

708 2(ax--6)/3—sin (2c) «sin (az)] 224. (33) 

(4) The Thomas—Fermi equation of atomic theory is, in dimensionless 
form, 

d 


CE at gr Watlie O(O}==l a id(co)==08) File.) (34) 


Referring to the general expression (6), it is seen that for this case 
—— —1/2 
Ay=D?, a=%, n=0, m=, ay=— x?! 


and all the remaining a, are zero. (d(x) in (6) is zero.) ; 

The generating function is clearly of the form (b-+-ax). If is put equal 
to 1 to satisfy the first boundary condition, the subsequent functions A, 
must be chosen so that they vanish when 2 is zero. Thus 


A,= & fecn .(1+ax)3/?.. da dx 
0/0 


=a-3/2 {2(1 +-27)1/2[z4/6 + 227/3 — 1/8] 
+(1+62?)/8 . log [z+(1+27)#7]} . . . . (35) 
| 1 (4 
bs Ber hon is cumbersome to handle in the evaluation of A,, and it 
is convenient to approximate by assuming that zis small. (35) may then 
be shown to yield the power series 

A, & a3? {428/34+225/5+327/70+...}. . . « « (36) 
This approximation may be used to extend the solution until ueiher 
terms are small compared with 2’. Proceeding in this way a power series 
is obtained for 4(x) 

Hla) = {1 +av-+42° 2/3 + 2a0%/5 
+a2/3+3a?x7?/70+2aa4/15 ...}. . ~ (37) 
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This series is identical with that obtained by Baker (1930) by a method 
of successive approximations. If the full form for A, given by (35) is 
retained, and an approximate form of A, is used a series will be found 
that converges more rapidly than (37). 

Tt is interesting to note that an alternative solution can be found by 
solving (34) for 4?” instead of ¢. In this case the parameters are 


dgee Le, a=5/2, n=—1, m= —5/3, d= — ee 


$7. CONCLUSION. 


It has been shown that it is possible to revert series involving certain 
fractional powers of the unknown variable, and that the coefficients in 
these series need not be algebraic. The application of reversion to 
non-linear differential equations has been found to be simple, and applicable 
to a wide variety of problems. It is particularly suited to problems with 
initial boundary conditions. It suffers from the presence of “ secular 
terms”? in common with other series methods (Minorsky 1947), and is 
not usually suitable when the non-linear terms are large. 


The author is grateful for permission to publish this work, from Standard 
Telecommunication Laboratories Ltd., and also from the British 
Thomson-Houston Company Ltd. where the work was commenced. 
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SUMMARY. 


The general equations of boundary layer flow along a curved surface 
S are set up in an appropriate system of orthogonal coordinates and the 
importance of the curvatures in the general case is made evident. 

It is shown that, for a general flow, the necessary and sufficient condition 
for the curvatures to disappear from the equations of flow is that the 
boundary layer should be along the surface of a cylinder or a plane. 

The form assumed by the equations when S is a surface of revolution 
is derived. 

In Part II. there will be a discussion of the flow in the vicinity of a. 
stagnation point on a general surface. 


THE COORDINATE SYSTEM. 
LET us denote the equation of the given surface S by 


MUS Clan U sy (i) Sele caret wa ge bd) 


This surface together with all parallelt surfaces constitutes a Lamé 
family of surfaces { 
(aed e2\==CONSt gate ee eo. F.. (2) 


Tt follows that two families of surfaces 
E(x, y, z)=const. and 7(%,y,z)=const. . . . . (3) 


exist such that the three families €, 7, ¢ form a triply orthogonal family. 
The intersections of the original surface with the members of € and 7 are 
the lines of curvature of S§ and the surfaces = const., y=const. are, 
fact, the developables formed by the normals along the lines of curvature 


on §]j. 


* Communicated by the Author. ; men? st 

+ See Weatherburn ‘ Differential Geometry ’, p. 158, “A surface which is 
at a constant distance along the normal from another surface S is said to be 
parallel toS”. 

+ Op. cit,, p. 207 or p. 218. Ex. 2. 

§ Weatherburn, op. cit., p. 211. Dupin’s Theorem. 

|| Op. cit., p. 207. 
“ 4 
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If h,, ho, hg are the length elements in our coordinate system then, 
by paneer h, is a function of £ only. We can, without loss of generality, 
take it to be unity, ¢ then representing distance measured normal to 8. 


THE Equations oF Morton. 


In terms of these coordinates, denoting by wu, v, w the corresponding 
velocity components, 


div V= 7 : marae w+ 5 (hw)+ se lahat) | 2.) 
Ow oO Lsfeo Ow 1 0 0 

curl V= > ey — 5 (hw) | iLae u)— ae a | get) gyno | 

5) 


The components of acceleration are 


du udu, v du Ou w oh, wv dh, , wu dh, 


Gi | hve | hy On) Olas hihgion.» buhe Olena ae 


BO th OP et OU ee es 6 
Of hy of hig 0n | OL’ “Tala On Raha oe ne OL re (6) 
aw u dw, vw, dw a8 ay _v* dh 

ot h, OE — hg oe hy Ot hs ca 

On making the usual boundary layer assumptions, with the additional 
assumption that h,, h, and all their first derivatives are O(1), we have the 
following equations of motion 


patie en 
and the equation of Sontintaty is 
rag | gent + Ou) | + SF =o, faite eid (0) 


The length elements h,, hy are not completely arbitrary ; the Lamé 
relations, in our case, lead to the following results :— 


(i. 1 oh, he I hy 
i. heen anc Ss = are each functions of € and 7 only, 


(ii.) A, and hy are each of the form F (€, n) +¢G(E, n), 


. 0 / Wh 8 /daae een eon . 
(ii1.) a2 (Fe tHE m) Ol ana —=(§) ‘ ° e ° * ® (11) 
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; Furthermore, if K, is the principal curvature on the surface €=const. 
in the direction of the parameter 7 (the curve ¢=const. on this surface) 


and Ky is similarly the principal curvature on the surface 7=const. in 
the direction of the parameter €. 


Wossl ron 


pee lirch, 1 oh, 
< hih, o€ 


id K,—=— — —— 
an : Teeiaal: 


(12) 


BI 
The surfaces €=const. and »=const. being developables, their other 
principal curvatures are zero, of course. 


Also if K, and K,, are the principal curvatures of the surface C—const. 
in the directions € and y increasing 


ihe Ce many a 
im h, ol’ Via | hee Sie acre 1 eis) 


Therefore, in terms of the curvatures K,, K,, the equations of 
motion (7), (8), and (10) may be written 


Cian UACUrme ty Ot. OU: - P 1 op. du 

Ot hy 0g Ponte oO et ae cae a (14) 

ov ud vd ov l op dv ; 

ET SEY: He Eade +K,u —K w= eon oCe (1 ) 
1 ou lav dw rs : 
ie By ss Bias (ian ke a Ss ae ME) 


where the values of h,, hy, K, and K, may, to the orders of approximation 
involved, be taken in the surface 8S. As in two dimensional theory, we 
must suppose p prescribed, in virtue of (9), by the flow outside and then 
equations (14)—(16) are the equations which serve to determine wu, v, w. 


UNIQUENESS. 


Considering the field of flow as a whole, the system of orthogonal 
coordinates we have set up is clearly not unique. For example, for an 
ellipsoid the confocal system is an obvious alternative. Nevertheless, 
since by Dupin’s theorem the curves of intersection of a triply orthogonal 
family are lines of curvature on each it follows that, however we set up 
the orthogonal coordinate system, the surfaces =const. and —const. 
must intersect the surface S in its lines of curvatures. Hence, no matter 
how we choose our triply orthogonal family to contain 8, the various 
possibilities for the families €=const. and 7—const. must all touch the 
families we have chosen along the lines of curvature of S. Since we are 
concerned in boundary layer theory only with the immediate vicinity 
of S our equations represent the unique form for an orthogonal coordinate , 
system, apart that is, from the choice of scales for € and 7. 
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The general flow along the surface of a right circular cone of semi-angle 
x makes an instructive illustration. The coordinate system we have 
developed consists of 


(i.) cones parallel to the given cone ; 
(ii.) cones orthogonal to the given cone ; 
(iii.) meridian planes. 


In cylindrical polar coordinates r, ¢, 2 we may express these families as 


é=¢, yn=zcosatrsina, f=zsina—rcosa, . . (17) 

so that h,=nsina—Ccosa=r, ha=hs=1, Aiins ee SMe 
ou 

Kr; K,=— 7, BIN." 5 eee 


Spherical polar coordinates R, 6, ¢ with corresponding velocity 
components U, V, W would be an obvious alternative system. In the 
comparing the two we should have 
u=W, v=U cos (6—«)—V sin (@—«), w=—Usin (@Q—«)+V cos (@—«)* 

(20) 
and at the surface of the cone 
Leo 1 0 0 0 Cae ig 
i, Devt Raine Ot ae I OR ea ee 
It can then be verified that the two sets of equations are identical. 


(21) 


DIscUSSION. 


The presence of the terms K, and K, in eqns. (14)—(16) makes manifest 
the importance of curvature in the general problem, though it should be 
noted that the significant curvatures are not the principal curvatures 
of S but the non-zero principal curvatures of the two orthogonal familes 
€ and y. In this respect, the theory differs from the two dimensional 
theory where curvature effects, so long as the curvatures are not very 
large or change very rapidly, disappear from the equations. We may, 
in consequence, in three dimensional boundary layer flow expect effects 
additional to those we might reasonably associate with the addition of 
an extra dimension. 

When K, and K, both vanish we have what is clearly an exceptional 
case in this respect and it is of some interest to determine the surfaces 
so defined. Although the vanishing of K, and K, refers only to values 
at the surface §, it can readily be established that they must in consequence 
vanish everywhere on the developables. Hence the developables 
€=const. and y=const. are planes. 

Furthermore, we may now establish that one of these families consists 
of parallel planes. For, suppose, if possible, that two members of 
€=const. are non-parallel planes. Any member of 7=const. must then, 
since it is orthogonal to both the chosen planes, meet them and must have 


* Since differentiations have to be carried out it is obviously insufficient to 
consider surface values of uw, v, w only. 
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a normal whose direction cosines are uniquely defined by the direction 
cosines of their normals. In other words, the family »=const. must 
consist of parallel planes. Hence we infer that at least one of the families 
é =const. and »=const. must consist of parallel planes. Suppose it 
1s y—=const. We can then without loss of generality express the equations 
of the two systems in the forms. 


aU | ng (ee DU oe ee Lee Leen ed 

Tae ROtey eee ene ae (3) 
Suppose the third family is expressed by 

LEN CORSO ea ot Avg tare SARA D AN de husk eed ala 2h) 


Then, considering the three surfaces with parameters €,, 71, ¢,, which pass 
through the point z,, y,, 2, we must have 


4=g, Ww=Af(S)+o(Es), Fey, %)=G, ~ « (25) 
and, since the surfaces are orthogonal 


(Ge aoe aa) ng a) 
O% Jame, _ \OY/yay, _ \O%emn |... (26) 


| F(E1) 0 


Hence the surfaces €=const. are, in general, cylinders with generators 
parallel to the z axis ; they include as a special case (when f(€)=0) planes. 

We may conclude therefore that it is only for boundary layers attached 
to cylinders and planes that the curvature effects vanish from the 
equations of motion. Moreover, in these cases, the boundary layer 
equations may be taken in Cartesian form. This is obvious for the case 
of the plane. For the cylinder it follows from eqn. (12) that h, is 
independent of 7 and h, is independent of € and therefore by condition 
(ii.) of p. 240 that in the surface S the length elements h,, h, can each without 
loss of generality be transformed into unity. Thus in the case of the 
cylinder the equations may be written in Cartesian form with € as distance 
measured around a right section, 7 as distance measured along the 
‘generators and ¢ as distance measured normal to 8. We thereby recover 
the standard result in two dimensional flow and obtain the result already 
justified in particular by previous investigators that the equations in this 
form are valid for three dimensional flow past a yawed cylinder. 

Finally we may obtain the forms taken by K, and K, when the 
surface S is one of revolution. The lines of curvature of 8 are then the 
meridian curves m and the “ parallels” p. Thus the surfaces ¢=const. 
are the meridian planes and hence 

Ky=0. 
The surfaces y=const. are the right circular cones generated by the 
normals to S along the curves p. Let us consider a particular surface 
n—=const. generated by the normals to § along a particular parallel »,. 
Then if R is the length of the normal to this cone intercepted between p, and 
the axis of rotation 


[K,|=1R. 
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SUMMARY. 


The results of the Bloch-Nordsieck method are derived in a covariant 
manner, based on a covariant form of the commutation-relations in 
momentum space. 


§ 1. 
In this paper we shall give a covariant derivation of the well known Bloch- 
Nordsieck (1937) results on the multiple production of low energy quanta. 
A full proof will be based on the interpretation of the field theoretical set 
of equations recently given by Yang and Feldman 1950. In §3 
a derivation of the same result from the Dyson-Feynman form of the 
S-matrix will be outlined. 

Though the results obtained here are not new, and have in fact been 
discussed by several authors, the importance of the Bloch-Nordsieck 
problem as the only problem in the theory of fields which admits an accurate’ 
solution seems to justify a general reformulation of the procedure. The 
simplification which enables one to find an accurate solution rests on the 
neglect of the recoil of the source particles, the motion of which is assumed 
to be given. The theory of fields contains a second problem of this kind in 
which the réles of Fermi- and Bose-particles are interchanged. Though 
this problem can be solved in principle, no manageable form of solution has 
- been found so far. If the solutions of both problems could be handled with 
equal ease they could be made the basis of a new approximation method. 
An approximation of this type would have the advantage that the trivial 
effects of higher powers of the interaction constant (unitarity of the S- 
matrix) would be considered in every step of the procedure, and only the 
genuine higher order effects expressing the intrinsic non-linearity of the 
system would be subjected to gradual approximation. It seems essential 
for such a programme that the Bloch-Nordsieck result should be given in a 
covariant form, though its immediate practical applications are of course 
non-relativistic. 


§ 2 
8 . 
We consider a system of coupled fields. The Bose field will be described 


by a wave operator ¢(x), and we shall restrict the first part of this discussion 
to scalar particles. The extension to more general Bose fields is 


* Communicated by the Authors. 
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straight-forward. The source particles are summarily described by a 


source density p(w); x denotes a point in space-time. The field equations 
are written in the Yang-Feldman form : 


P(x) =$"(a)+ fda" p(a’)D,(w—a') 
=$N(a)+ fda! p(x’)D,(a—2'), 


where the functions D,(«) and D,(x) are the retarded and advanced solutions 
of the Poisson equation 


ee) (y= O(a) sea weir arene rrih(7) 


{ . 

Since both 4 and ¢°“t are solutions of the homogeneous wave equation 
(L?—?)¢"=—0 and, therefore, satisfy identical commutation relations of 
the form 


peo (es l= Dewy. aN, LO BT (2) 
they must be connected by a unitary transformation 
PM (ae) = P(x) + f(D,—Dg)p(2' da’ =S 1638, . . . . (3) 


and it has been shown by Yang and Feldman that the unitary matrix § 
thus defined is identical with the Feynman-Dyson S-matrix. 

Since ¢ is a c-number—the motion of the source particles is fixed—d 
differs from ¢" only by a multiple of the unit matrix. D,—D, is an odd 
invariant solution of the homogeneous wave equation, and actually 


D,(x)—D,(x)=D(2). 


The further calculation is considerably facilitated by transforming to 
k-space. This gives 


Cut) — ap | dies 2 —peygon (ket 
= yp | Cen? NEBME)—pl bel etd) 


where the first equation is a definition of the Fourier transform h(k). 
Use has been made of the Fourier representation 


Die ae | dk3(k2—p22)e(kethe 
of the D-function; <(k) is --1 according to k)20 and p(k) is defined by 


th af ; 
ree. we AS 
pl”) = oi | akp (ye (5) 
In the following the abbreviation 34(4)= —p(k)e(k) will be used. 
The energy momentum 4-vector for the incoming wave is given by 


‘ 1 ; inf J.\+in/] 
pues tae | dk (k2—p2)k,b™(k)+ 6h), 
Fred et Man 
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and a similar expression holds for the outgoing wave. The number of 
quanta in the incoming wave is 


Nina Re dd (k2— 2) (ke) +Gi8(ke) es ( 
ko> 


(27)? 
In the present case we want to consider an experiment in which quanta are 
only observed over a limited range of frequency. We shall therefore 
introduce an interval 4 in k-space which contains only this range of fre- 
quencies. It is then easily seen that N,—resulting from equation (6) 
by integrating only over the interval 4—is still an integral positive number. 
This follows from the commutation relations 


8(k’2—p2)3(k2—p2) [b+ (he) bk’ )] = —(2m)93(k— ki )e(k)B(k">—p2) (7) 


for the free fields, which can be derived from equation (2). Using (7), one 
obtains 


5(k?—p*)[N ah(k)] = —9(k)0(k?—p?) 


8(k? — w?)[N nh*(k)]= + 6*(k)8(K? —p?) 
for any ke A. i 
Since we want to determine the probability for the production of a 
certain number » of quanta in 4, we shall assume that the incoming wave 
does not carry any quanta at all in 4. Then if yp is the eigenvector 
(| Xo |?=1) describing such a state, we have 


O(k?— ph) vos" coc cs ae te 
for every ke 4. For the outgoing wave y;, is the corresponding vacuum 
state and the basis-vectors for the 1, 2, .. . n particle states are 

xi (hy “wr | phy vo(ki—p?), 
x2(hy ko) ai Pe (ks )pr(heo)* x3 (hA7—p?)8(k3—p?), 
Xn(hy - . . ky)= Bie a +2 PN len) xaLTB (ki). 


The most general n-particle state can then be described by a Hilbert 
vector, 

e: Wes 

Fi= gap | , 2 eH he Ober (ey. . ky xglky ++ ky)e . . (9) 


where wi(k,,...k,) is a symmetric ‘ Schrédinger-function ” and the 
index r indicates one particular distribution of the quanta over the element 
A of k-space. To normalize F” such that F7N,F”—n we have to make 


7 
( 277) 8h 


| deren] dk, L15(?—p?) | w, P= 
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The proof of this is analogous to that for the well known case of oscillator 
eigenfunctions and can be carried out simply by using the method developed 
by Becker and Leibfried (1946). We shall use a complete set of orthogonal 
‘functions w’,, i.e. such a set that every symmetric function f,(k,,...k,) 
defined in 4 can be expanded in terms of the u,. The w’s will then 
satisfy the completeness relation 
1 1 
ene = pie oe dy (10) 

The symmetrized 5-function on the right-hand side of this equation can be 
replaced by the ordinary 5-function whenever it occurs multiplied by a 
function symmetric in the k,,...k,,. 

The probability amplitude for the creation of n particles in a state r is 
now given by 


ETTS(he—p2)uh*(hey . «I, )eah (del... 


F vente: 


1 1 . - 
Oil et oe 7 C p22 
(Fixe) = Gyaym X api) , Mea + |, eat TOUE—H2) 


XU (hy - « - Ky)(Xo( Pes) +86(hy)) - « - (Phy) +86(h,)) Xo): 
On account of equation (8) only the product term 5¢4(k,)d4(kg) . . . dA(k,,) 
gives a non-vanishing contribution. Since we are not interested in the 
particular distribution of the quanta inside the interval 4 we have to sum 
the squared amplitude over all the distributions defined by wu’. Thus, by 
using the completeness relation (10), 


: 1 mir , : 
“EF iXxo We obemae PYerel On Aye cee (11) 
, FS 1 ‘ i 2__ 2) 1 §h(k) |2 pay fee 12 
where = ae | ahah 12) (OD) tae (12) 


Since the total probability must be one we have 


Re Te ae vs ye (13) 
n is the average number of quanta produced: n=2nw,(4). It can also be 


nr 

found by considering the vacuum expectation value yyN°"'yo of N a 

Since the interval A is arbitrary, this—in a sense—is a generalization of 
‘the Bloch-Nordsieck result. It shows how redundant information, e.g. 
the distribution of energy and momentum over the quanta in A, can be 
eliminated by the use of the completeness relation. Another link with the 
case of infinitesimal intervals considered by Bloch and Nordsieck lies in the 
addition theorem of Poisson distributions : Let the interval 4 be divided 
into two subintervals, 4, and 4,. Then n=7,+7, and 

(y+Me)" 


w,(4)= nd ae exp{ a (n, +7») }= = On (A oye 2). 


The validity of the Poisson distribution in two subintervals therefore 
insures the validity of a Poisson distribution for the whole interval, and the 
proof can be easily extended to cover the case of an arbitrary number of 
subintervals. 
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It remains to define the source function p(k) for the particular case where 
a source particle is scattered at ~=0. Transcribing the above argument to 
a Maxwell field, we replace p(w) by a vector source function j,,(%) of a point 
electron, viz., 


jee] pul2)8(e—7P(2) dr 


where, e. g., p,(7)=p, for 7<0 and p,(7)= =Die for 7>0. The sudden change 
in momentum imposes the restriction that in applying the result to a real 
scattering process one should limit the admissible frequency range to 
frequencies <1/r7, where 7 is the effective time of collision. Otherwise a 
sudden change of momentum would be too rough an approximation. 
According to equation (5) we have for the Fourier transform of j,,(x) 


io=«( G5 wa) 


The average number of quanta can be determined from equation (12) : 


ee | aks(i2—p2) 
* Grd n ‘ 


(pe) _ (p’e) 
(pk) (pk) 


where e is a polarization vector. 
In the case of mesons p(a) will be of the form 


ple)=g | o(7)8(e—rp(r))dr, 


where o(7) is a function describing an internal degree of freedom of the 
source particles, as e.g. their spin or isotopic spin. Both pando can change 
at 7=0. In this case the Fourier transform of the source density is given 


by 
er at=)_ ott} 
Ali = Tm Dh J" 


n can be easily determined from equation (12). 


We shall now briefly indicate the derivation of the above results with 
Dyson’s S-matrix 


S=1—i| pt)p(u)+! 2)p(2)-+...  . (14) 


Here $(1) means 4(x) at the point 21, P is the operator arranging the factors 
in the product in or ie of time, and the integrations are to be extended over 
the whole x,,...2, space. The problem consists in finding the matrix 
elements of S which lead from the vacuum state to a state where n particles 
are present. To this matrix element we shall get contributions from the 
n+ Ist, n+ 8rd, ... terms in the expansion (14). A typical term of n+ 2/th 
order will contain n factors 4(1) which create the final particles, J factors 
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which create / virtual particles, which are finally annihilated by the 
remaining / factors. In calculating the matrix element one has to keep 
in mind that one gets a factor 
n+ 21 
n 


corresponding to the number of possibilities of picking n 4’s out of the total 
n+2l ¢’s. A further factor n! arises from the permutations of the n 
created particles among the n creation operators (1). The matrix 
element will therefore be of the form 


O[S[m)=(—a)"fP{CO|b] ky)... OD], de(1) . . « p(n) 
x | 14+ ae (0 | d(n+1)d(n+2)|0)4 se ih Mee ee (1) 


where Dirac’s notation for the matrix elements has been used. It should 
be noted that the expression in the square brackets is independent of n. 
The relative probabilities for the production of a given number of particles 
are thus given by the first non-vanishing approximation. Integrating the 
squared matrix element over all energy momentum vectors of the quanta, 
we obtain 


1 1 
KOSI") P= am Xa], ge Modeled» Lok 


x 118(k?—p2) [ TP. 


The 1/n! arises from the fact that one wants to integrate all the k’s over the 
same interval. Since the particles are not distinguishable, one has to 
divide by n!. [ ]?, which is the square of the bracket in equation (15), can be 
determined as before with the aid of the unitarity of the S-matrix. It 
gives only a normalization of the probability. The result is again 
Leal by 
[Ko] S|ny[2= Se", 


n! 


‘with n, as before, defined by equation (12). 
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ABSTRACT. 


Measurements with a penetrating shower set at sea level showed an 
apparent mean free path for absorption of the secondaries of penetrating 
showers of 405-430 gr./em.? for lead and 350-+-40 gr./em.? for iron. 


INTRODUCTION. 


MEASUREMENTS for the absorption of penetrating shower secondaries have 
been made previously for lead for thicknesses of 100 gm./em.?— 200 gm./em.? 
(Walker 1950). A comparison of the differential absorption for lead and 
iron under a large thickness of lead (150 gm./em.?) has been made by 
(Piecioni 1950). 

The measurements reported below give results for the absorption of the 
secondaries in the range 60 gm./cm.?—450 gm./cm.? for lead and 60 gm./cm.* 
—260 gm./cem.? for iron. 

These investigations were started at a time when it was considered that 
the penetrating particles of a shower were predominantly 7 mesons. 
Later results have shown that the particles capable of transversing 20 
centimetres of lead will be composed almost equally of protons and 7 
mesons (Camerini et a/. 1950). Thus such an investigation does not throw 
direct light on the 7 mesons interaction. 

However, the background data provided by such an investigation of 
the effect of lead absorber thickness on the counting rate of a ‘“‘ penetrating 
shower” set and the further data obtained on the penetration of the 
shower secondaries of an ‘‘ extensive shower ’’, appeared to make the 
publication of a brief report of my investigation worth while. 


: EXPERIMENTAL. 

Measurements were made on penetrating showers at sea level with the 
apparatus shown in fig. 1. 

The top six counters (marked either A, B, or C) and the lower ten counters 
marked either D or E were of active area 3-5 em. 50cm. Each counter 
in the top tray was separated from its neighbour by 2:5 em. of lead. The 
tray was covered by 15 em. of lead, the side shielding consisted of 10 em. 
of lead, and there was a thickness of 2:5 em. of lead immediately beneath 
the tray. 


* Communicated by Professor C. B. A. McCusker. 
+ Now at Department of Physics, University of Rochester, 
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The absorber 8, over the lower ten counters was variable. The side 
shielding again consisted of 10 cm. of lead. 

Coincidences (A B C D E) were recorded. In order to record penetrating 
showers simultaneously accompanied by air showers, the pulse X; was 
taken from coincidences between three extended trays, and coincidences 
ABC X;, A BCD E X, were recorded throughout the experiment. (The 
three extended trays each of area 750 cm.2 were simultaneously used in 
another experiment, in which they were covered with thicknesses of lead, 
water, carbon, iron, up to one cascade unit, thick.) 


Bios! 


oe 


SS ERO 
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RESULTS. 


Fig. 2. plots the logarithm of the rate A BC D E against arraise Ms 
the absorber S, in g./em.2 of lead. It will be seen that the results can : ‘ 
represented by an exponential absorption for lead with mean free e 
405+30 gm./cm.2. (The results were corrected for barometric variations : 
a pressure of 1015 millibars using a barometer coefficient of —1 — 
per millibar.) A measurement was made with an absorber thickness 5, 


Counts ABCDE per hr. 
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composed of 7:5 cm. of lead directly over the counters, and 48 gm. /om.? of 
carbon above the lead. The observed rate was 0-870-05 per hr. 
Comparison with the lead absorption curve shows that this rate corresponds 
to a lead absorber thickness of 90+40 gm./cm.?. 

Measurements made with absorber 8, composed of 260 gm./cm.? of iron 
gave an observed rate corrected to 1015 mbs atmospheric pressure, of 
0-52-03 ¢.p.h.r. corresponding to a lead absorber thickness of 300 
gm./em.2 thus giving the stopping power of iron relative to lead 
as 1-15-L0-15. 


0 
0 aaeatah00 200 300 400 500 600 


Absorber S, in gms./em.? 

Fig. 2 shows the counting rate A BC D E plotted against absorber thickness 
in gm./em.?, The triangular point at 260 gm./cm.? refers to iron absorber. 
The full line is the best curve for lead absorber, and the dotted line is the 
best for iron absorber. 


Table I. gives the rate of accompanied showers. As the absorber over 
the extended trays was varied during the course of the experiment we 
compare the ratios 

ABCDEX, 
ABCX; 

Including all results for small absorber thicknesses we obtained a value 

of 
ABCDEX, 
ABCX, 


of 0:71-+0-10. Ifthe absorption of the accompanied showers was the same 
as for non-accompanied showers we should expect a ratio of 0-31-40-04 
for 450 gm./cm.? of lead absorber. The observed value of 0-55--0-10 gives 
an indication that the accompanied penetrating showers are more 
penetrating than the non-accompanied, though the statistical weight is poor, 
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Discussion and Comparison with other Experiments. 


Similar penetrating shower measurements have been made by’ Piccioni 
(1950), Walker (1950). The number of penetrating secondaries, however, 
necessary to cause a count varied in the two experiments. The data of 
Walker (1950), also Piccioni (1950) show that the repartition of shower size 
under various thicknesses of absorber does not appear to alter. Accord- 
ingly the absorption curves obtained in the different experiments should be 
similar. 

The values obtained for the apparent mean free paths in Lead by 
Piccioni and Walker, are respectively, 360+40 gm./em.2 and 405+110 
gm./em.”. Our value of 405+.30 gm./cm.? isin good agreement. Further 
Piccioni obtained a ratio for the stopping powers of Iron and Lead of 
1-:02+0-13 in agreement with our value of 1-15--0-15. 


TABLE I, 


Data on Penetrating Showers Accompanied by Extensive Air Showers. 


Thickness of S Number of Number of 
Absorber Time (hrs.) Counts Counts 
A p.G x, ABCDEX, 

26 gm./cm.? of lead 97 : 3 3 

78 gm./cm.? of lead 140 9 vk 

105 gm./cm.? of lead 128 11 5 
78 gm./cm.? of lead 
plus 50 gm./cm.? of 

carbon 339 18 14 

210 gm./em.? of lead 284 21 12 

315 gm./cem.? of lead 137 8 3 

450 gm./cm.? of lead 500 31 ii 

260 gm./em.? of iron 444 29 18 


The interpretation of these results is not straightforward, but it is now 
clear from the results with both cloud chambers and photographic plates, 
(Brown and Mckay 1950, Camerini e¢ al. 1950) that the high energy 
secondaries of penetrating showers are a mixture of protons, neutrons and 
' 7-mesons. 

At any level in the absorber there will be simultaneous absorption and 
production of ionizing particles, the measured absorption thus being only 
indirectly related to the actual absorption processes taking place in the 
material. 

Stress has been laid (Piccioni 1950, cf. Greisen 1950 for discussion of 
Piccioni’s results) on the fact that if nuclear processes predominate for 
absorption and production the stopping powers of iron and lead should be 
in the ratios of the geometric cross sections 7. e. 1-55: 1, and not in the 
ratios of the stopping powers for ionization losses of 1-1 : I. This, however, 


aa) 
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is too simple an approach to the problem. In the presence of some degree 
of transparency such as would be expected for nucleons with energies of the 
order of 1 Bev, and for 7-mesons in the 100 Mev—1 Bev range, the ratio of 
the cross sections for nuclear interactions would be considerably lower. 
In addition ionization loss processes will play some role and a ratio of the 
order of 1-2—1-4 would not therefore be in contradiction with present views 
as to the relatively strong interactions of 7-mesons and nucleons. 

It is to be expected that the nucleons forming part of the extensive 
showers will have a harder energy spectrum than those nucleons not accom- 
panied by extensive showers, and thus those nucleons should give rise to 
higher energy events with greater penetrating power. 
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SUMMARY. 


Schlieren photographs have been taken to show the supersonic flow past 
a sharp wedge placed at an incidence large enough to give an Fan 
region on one side. In addition to the flow phenomena predicted b 
inviscid theory, the photographs show a weak shock wave Selbainteiy 
upstream of the expansion region. Various possible explanations of this 
are considered and it is shown that the shock wave is caused by the Shan) 
bluntness of the leading edge, although the thickness of the leading edge is 
only about 8. It cannot be explained by considering the growth of 
the boundary layer, neither is there any evidence of boundary-layer 
separation. 


§ 1. INTRODUCTION. 


AccorDING to inviscid theory, the supersonic flow past a perfectly sharp 
wedge consists of two independent regions separated by the wedge. If 
the angle of incidence of the wedge is sufficiently large there is a centred 
expansion on one side and a shock wave on the other. it was expected 
that this theoretical flow pattern might not be found in practice, either 
because of viscous effects or because of bluntness of the leading edge. 
Some experiments were therefore made to investigate the supersonic flow 
past a wedge at a high incidence. 


§ 2. NoTaTION. 


M, Mach number. 
p, pressure before shock wave. 
Ap, rise of pressure in shock wave. 
R, Reynolds number. 
s= Ap/p, shock-wave strength. 
T, absolute temperature. 
t, thickness of body near leading edge. 
u,, velocity outside boundary layer. 
x, distance from leading edge. 


* Communicated by W. A. Mair, 
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y, distance from axis. 
8,, displacement thickness of boundary layer in compressible flow. 
5*, displacement thickness of boundary layer in incompressible 
flow. 
y, ratio of specific heats. 
6*, momentum thickness of boundary layer in incompressible flow. 
uu, coefficient of viscosity. 
v, coefficient of kinematic viscosity. 
o, Prandtl number. 


§ 3. EXPERIMENTAL DETAILS. 


The intermittent supersonic wind tunnel, the two-mirror schlieren 
system, and the technique for measurement of humidity have already been 
described by Bardsley and Mair (1951). The mean Mach number in the 
working section during the present experiments was 1-965. An alter- 
native optical system which produced a magnified image was used for 
figs. 6-8 (Pls. XI. & XII). In this system the object lens was removed 
and the second mirror arranged to give an enlarged image of a selected area 
of the working section on the screen. Measurements of the humidity of 
the air passing through the tunnel showed that the water-vapour content 
was always less than 1 part in 2000 by weight. 

The wedge used was 1-27 cm. long in the direction of the stream, had an 
apex angle of 10-0° and extended across the whole width of the tunnel. 
It was made of high-carbon steel, ground and lapped to make the leading 
edge as nearly perfect as possible, with a constant apex angle. An estimate 
of the thickness of the leading edge was obtained using a shadowgraph 
technique. The profile of the wedge was illuminated and a short-focus 
convex lens was used to form an image magnified 25 times on a screen. 
Several typical parts of this image were photographed before and after the 
experiments. Before the experiments about 90 per cent of the span had a 
profile similar to that shown in fig. 1 (Pl. IX.) ; the remaining 10 per cent 
was similar to that shown in fig. 2 (PI. 1X.). After the experiments about 
50 per cent of the span had deteriorated to the condition shown in fig. 2, 
the remainder being unaltered. The deterioration was probably caused 
by dust in the air stream, although the total running time of the tunnel 
was only about 2 minutes. 

Fig. 9 is a diagrammatic representation of the indentations shown in 
figs. | and 2. The average depth of the indentations (h) was measured on 
the photographs. Assuming that the wedge angle was constant and that 
the thickness at a point A was very small (2 4), the thickness at a point B 
was calculated. This thickness was found to be 5 u for the parts of the 
edge in better condition (fig. 1) and 12 p for the deteriorated parts (fig. 2). 
The same method and assumptions were used previously on a similar wedge 
and gave good agreement with the thickness as measured by a calibrated 
microscope, 
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The wedge was supported in the wind tunnel by two arms on the down- 
stream side. The position of these arms was altered, as required, to change 
the angle of incidence of the wedge. 

The knife edges in the schlieren system were horizontal and were 
arranged so that in all the photographs increasing density in an upward 
direction was represented by a dark region, 

Preliminary visual observations, using a mercury-vapour lamp, showed 
that the flow on the lower side of the wedge was rather unsteady when the 
deflection angle at the upper surface was greater than about 7°. This 
unsteadiness appeared to be due to the interaction between the shock wave 
on the lower side of the wedge and the boundary layer on the bottom wall 
of the tunnel. Because of this unsteadiness a short-duration flash 
(effective exposure about 5 microseconds) was used as the light source for 
photography. Photographs of the flow are shown in figs. 3-5 (Pls. X. & XI.) 
for positions of the wedge giving stream deflection angles (on the expansion 


Fig. 9. 


=e 


Indentations in leading edge. 


side) of 0,5,and 10°. (‘The possible error in these angles is on) At the 
Mach number of these experiments (1-965) the maximum deflection at an 
attached shock wave is 22-3°, so that with this wedge (10° apex angle) the 
deflection at the expansion could not exceed 12:3". 


§ 4. DiscUSSION OF THE PHOTOGRAPHS. 


The arms supporting the wedge continue downstream beyond the field of 
view ; the position of the back of the wedge is revealed by the expansion 
region on the lower side. The wake is clearly visible in figs. 4 and 5, and 
shock waves can be seen where the direction of flow changes to that of the 
wake. The bump that appears about two-thirds of the way along the 
upper surface of the wedge is a small piece of rubber sealing-material goons 
of the glass side plates. In figs. 3-5 the region representing the shock a 
on the high-pressure side of the wedge is fan-shaped in oe a a 
is caused by curvature of the shock wave as it enters the turbulent 
boundary layers on the glass side-walls of the tunnel. 
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In each photograph the flow past the wedge is similar to that predicted 
by inviscid theory, but in addition there is a shock wave immediately 
upstream of the expansion region. This must be caused either by viscosity 
or by the slight bluntness of the leading edge. Three alternative explan- 
ations are considered below. 

Firstly, consider the rate of growth of the boundary layer. In the 
experiments the shock wave occurs for a stream deflection angle (measured 
at the surface of the wedge) as large as 10°. The boundary layer is 
therefore considered at the point where its rate of growth is equal to 10°. 
Blasius’ solution (1908) of the laminar boundary-layer equations for 
incompressible flow past a flat plate shows that the displacement thickness 
is 


6% =1-72(va/u,), 
and the momentum thickness 6* is given by 
6* /5* =0:386, 


where v is the kinematic viscosity, wu, is the velocity outside the boundary 
layer and «x is the distance from the leading edge of the plate. It follows 
from Howarth’s solution (1948) for a flat plate (assuming o=1 and pw « T) 
that for compressible flow the displacement thickness is 


8, =5¥ Eee (145) |. 


Hence in our case an approximate value of the displacement thickness of 
the boundary layer on the surface of the wedge is given by 


8, =4(va/u,)}. 


The rate of growth of the boundary layer is then 


$ 
Zu n2(25) —2R-3, 
dx Ux 


where R=w,2/v, the Reynolds number at the station a. 
At the point where the rate of growth is equal to 10° (0-2 radians), 
ay —0-2=2R-4, 
du 
therefore R=100. 

The boundary-layer theory is not strictly applicable for such a small 
Reynolds number. However, since the theory would make d8,/da— 0 
as «> 0 it follows that the actual rate of growth near the leading edge 
must be smaller than that given by the theory. Hence at the point where 
the rate of growth of the boundary layer is 10°, 2 is not greater than the 
value given by R=100. The stagnation conditions are atmospheric, so 
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that the Reynolds number in the test section can be expressed as 
1-34 10° percm. Therefore when R=100, 
<p UU eine 
EET Sx 10=* cm... 
and at this point 
6,=4R 4a 3 3x 10-4 em. 
Thus if there were a point at which the boundary layer grew so rapidly 
that the sign of the stream deflection was consistent with the presence of a 
shock wave, the values of x and 8, at this point would be of the same order 


Fig. 10. 


Flow with separation of boundary layer. 


of magnitude as the thickness of the leading edge itself. This suggests that 
the effect of the boundary layer on the formation of the shock wave in front 
of the expansion is small compared with the effect of the bluntness of the 
leading edge. 

Secondly, consider the possibility of boundary-layer separation. This 
is a well known phenomenon in the case of flow past blunt-nosed bodies and 
has been observed, for example, by Liepmann (1950) and by Holder, 
Tomlinson and Rogers (1949).. In their experiments the bodies were at 
zero incidence and it was found that for sufficiently blunt noses there was a 
separation of the boundary layer near the nose, followed by re-attachment 
with an oblique shock wave at the point of re-attachment. If this type of 
flow occurred on the expansion side in the present experiments the pattern 
observed would be similar to that shown in fig. 10. In particular, a shock 
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Fig. 1. 
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See rige?: : 
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Fig. 5 
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Stream deflection angle 10°. 
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wave would be visible on the downstream side of the expansion. No 
shock wave can be seen in this position in figs. 3-5, but to investigate this 
further the schlieren system was modified, as described earlier, to give an 
enlarged image. Photographs produced by this system are shown in 
figs. 6-8 ; the scale of these is such that the wedge chord is approximately 
three-quarters of the diameter of the field of view. These show the flow 
for stream deflection angles of 0, 5, and 10°, and correspond to figs. 3, 4, 
and 5 respectively. Again, there is no sign of a second shock wave that 
would indicate boundary-layer separation. Moreover, a theoretical 
investigation by Mair (not yet published) has shown that a shock wave 
having a strength (s=Ap/p) as small as about 0-01 should be detectable 
with this schlieren apparatus. The deflection of the stream at a shock 
wave of this strength is only about 0-2°. 

Lastly, consider the effect of the finite thickness of the leading edge. 
The inclination of the shock wave at a large distance from the leading edge 
was measured and found to be about 31-7° in all the photographs. This 
corresponds to a shock-wave strength (s) of 0-08. Since the possible error 
in measuring the shock-wave angle is about +4°, the corresponding error 
in s is +0-035, so that only the order of magnitude of the shock-wave 
strength can be found by this method. As already mentioned, a shock 
wave having a strength of this order should be easily visible in the 
photographs. 

Now consider the theoretical relation between the shock-wave strength 
and the thickness of the leading edge of the wedge. For inviscid flow, 
Friedrichs (1948) has shown that the strength of the shock wave from the 
leading edge of a thin two-dimensional body with a rounded nose is 
inversely proportional to the square root of the distance from the axis, for 
large values of this distance. This can be written in the dimensionless 
form 

s=k(t/y)?, 

where k is a constant (for geometrically similar bodies), y is the distance 
from the axis, and ¢ is a length defining the size of the body, taken as the 
thickness near the leading edge. Applying this result to the photographs 
of the flow past a flat plate with a blunt nose, obtained by Holder, 
Tomlinson and Rogers (1949), the constant k is found (by measuring the 
shock-wave angles) to be about 3. If it is assumed that the irregularly 
shaped edge of the wedge (figs. 1 and 2) gives a flow at large distances similar 
to that produced by a two-dimensional blunt nose of thickness ¢, the 
equation s=3(t/y)' can be applied to the present experiments. 

Putting t=8x10-4cm. (an average value) and y=6cm., we get 
s—0-035. Since the value of y at the upper boundary of figs. 3—5 is only 
about 5 or 6cm., this calculation shows that a shock wave having a 
strength of the same order of magnitude as that observed can be explained 
by the known slight bluntness of the leading edge, without any consider- 


ation of viscosity. 
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§ 5. CONCLUSIONS. 


It has been shown that the difference between the observed flow past 
a wedge (as sharp as possible) and that predicted by the idealized theory 
can be explained by the bluntness of the leading edge. There is evidence 
that viscosity has no significant effect on the formation of the shock wave 
upstream of the expansion region. 

If experiments of this kind were made with a leading edge of greater 
thickness, the strength of the shock wave for a given value of y would be 
greater. The effect of the growth of the boundary layer would then be 
even less important, and only if the boundary layer separated could 
viscosity have any significant effect. 

It may be noted that in these experiments the Reynolds number, based 
on the thickness of the leading edge, was very small compared with that in 
most other experiments on flow past blunt-nosed bodies. It is suggested 
that this difference in Reynolds number is important and may explain the 
absence of boundary-layer separation in the present experiments. 
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XXIX. Matrix Elements for Octupole Radiative Transitions. 


By A. P. Frencu, 
Cavendish Laboratory, Cambridge *. 


[Received December 11, 1950.] 


DuRING the analysis of some recent experiments on the angular correlation 
between alpha particles and gamma rays (Barnes, French and Devons 
1950) it was necessary to find values of the matrix elements for certain 
octupole transitions. These were derivable from general formule 
obtained by Fierz (1949) by group-theoretical methods. Although Fierz 
deals with multipole radiation of any order, his formule do not cover 
all possible transitions. If j represents the total angular momentum of 
the initial state, which goes over to a state of angular momentum j+ 4j 
with the emission of radiation of multipole order 2/, then the cases treated 
by Fierz are limited to 4j=+/ or +(J—1). It was thought to be of 
interest to extend calculations to all possible values of 4) for the particular 
case of octupole radiation, since radiation of this order (or even higher) 
seems to appear rather frequently in nuclear problems. 

The method used was a simple extension of that applied by Condon 
and Shortley (1935) to the problem of quadrupole transitions. We assume 
that the initial state is characterized by j, by its component m along some 
specified axis (which we shall take as the polar axis Oz in spherical polar 
coordinates, and by other observables «, with which we are not imme- 
diately concerned. The final state is characterized in a similar way by 
(j',m’, «’). The matrix element of electric octupole moment for a 
transition between the states may then be written as follows : 


(ajm | rrr | «’j’m’) 
— oy 3B (ajo | x | 0079’) (009m! | ox | el9"' om”) 


Coe fat rad aaa 5 21 
~~ (009m "”” | ¢ | o.'9’m'). 
For given values of («’’, j’’) and («’’’,j’”’) this summation has for all (m", m’"") 
a common factor G(j’’, j’”’) given (in the Condon and Shortley notation) by 
G(9”, He Ve= (09 3 / gy he 4a : / : gee ae yo : / E a9). 


In order to perform a final summation over the (j’’, j’”’) and so obtain the 
octupole matrix element in manageable form, it is necessary to relate the 
various G(j’’, 7’). This can be done through the requirement that rXr 


shall be zero. If we put 
r=ir+jy+kz, 
R=a—1y, 
then one must have zR—Rz=0. 


* Communicated by Professor D, R. Hartree. 
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This equation can be multiplied through from the left or from the right 
with z or R, and one can then demand that matrix elements such as 
(ajm | 2z2R—zRz | «'j’m') shall be zero. The matrix elements of R and z 
in the (j, m) scheme are tabulated (Condon and Shortley, p. ae so it is an 
easy matter to obtain a set of equations connecting the G(j”, 7”). Once 
this is done, a matrix element of octupole moment can be expressed in the 
form 
(ajm | rrr | a'j’m')=g(j, m, 7’, m’) XG, 

where G is the same for all transitions m— m’ having («, j) and («’, 7’) the 
same. g is then a definite function of the j’s and m’s and can be found. 

The radiation field is obtained from the electric moment by writing 

E(6, 6) =const X (8499+ Popo) . (%jm | rrr | «'7’m’) XoXo, 

where rp, 9, pp are the unit vectors of a spherical polar system. The 
double dot product here denotes a scalar product of ry with each of the two 
right-hand members of every triad of vectors appearing in the octupole 
moment. The form of the field depends, of course, only on 4m. The 


possibilities are listed below for convenience, normalized to an outgoing 
intensity of 47/7. 


E(0) = suey sin A(5 cos? 6— 1) — V"1(0) say, 


4 
1 ; 1 
E(+1)= <= 3 (8 cos 6(15 cos? @— 11) -epo(5 cos?6— 1) je=” = cae » 


» Sin 6(3 cos? 6B—1) +i, sin 26Jet*” = 


: to ~ 
E(+3)= vi (8, sin? 6 cos 6 Lig, sin® 6Je=2"” = L(+3) ,, 


The matrix elements for all possible types of electric aan transition 
are set out below. The corresponding formulae for magnetic transitions 
are obtained by rotating E through 90° about ry in the usual way. The 
matrices are, of course, Hermitian, but for convenience in use a complete 
tabulation is given. The coefficients A, B, C, etc., are the quantities G 
mentioned above. 

Aj=+3 
(J, m[E|j+3, m+3)=F fA {(jtm+1)(j+m4+2)(j+m+3)(j-+m+4) 

; x (jm-+5)(j m+ 6) }L( 3), 
(j,m|E|j+3, m42)=tAV/{(j+m+1)\(j j—m-+-1)(j +m-+2)(j7 --m-+3) 
x (gm-+4)(j +m-5)}L( +2), 
(j,m | E|j+3, m4£1)=F JAV/((j-+-m+1)(j—m-+ 1)(j-+-m-+2)(j—m+2) 
Xx (ga-m+3)(j+-m+4)}L(+1), 
(j, m|E|j+3, m) =—AvVt((j+m+1)(j—m+1)(j-+m+2)(j7—m-+2) 
x (J-+m-+3)(g—mM+3)}L(0). 


Aj 8 
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(j, m|E|j+2, m+3)—3B V (JF m)(j m+ 1)(j Lm+2)(j -m-+3) 


x (j4m+4)(j +m+5)}L(+43), 


(j,m|E|j+2, m+2)== $B(2jF 3m) /{(j--m+1)(j+m+2)(j--m-+3) 


xX (j-m-+4)}L(+2), 
(7, m|E|j+2, m41)=}B(j—3m)V/{(j+m+1)(j—m+1)(j--m+2) 
X (gm+3)L(+1), 
(j,m|E|j+2,m) =—3Bmy/{(j+m+1)(j—m+1)(j+m-+2) 
x (j—m+2)}L(0) 
Ai 


(j, m| E]j+1, m43)= +40+/{(jFm—1)(jFm)(j+m+1)(j+m4+2) 


XK (jy -m-+3)(j -m+-4) }L( +3), 


(j,m|E|j+1, m42)=— 3C(jF 38m—2)/{(j=m)\(jLm+1)\(j+-m+2) 


Xx (g+m-+3)}L(+2), 
m(27—1) 


ye | 
(j, m|E|j+1, m+1)= a 5 me? | Vigim+1) 


x (j--m-+ 2)}L( +1), 


(j,m|E|[j+1,m) =— of ee ADA) m | VG m+1)(9-+m-+1)}L(0). 


Aj=0 


(j, m|E|j, m43)=—}DV/{(j-m—2)(jFm—1)\(jFm)\(j+m+1) 


x (jm-+2)(j m+ 3)}L( +3), 


(j, m|E[j, m£2)= FAD Lim) V/{(GFm—l( jm j+m+I(j+m+2)} 


xL(+2), 
eee eae etre 
(J, mJ Ej, mae1)—4D | mm | VGFmG m+) 
x L(+), 
(2487-1) a a 
ip 


(j, m| E|j—1, m43)=F EV {(Gbm+VGem+ AIF mMGFmM—1) 


) 
x (j-Fm—2)(jFm—3)}L( +3), 


m|E|j—1,m+2)=—4 B(j£3m+3)v/ ((gtbm+1)\(jFm)(jFm—1) 


(); m|E|j—], m+l)=+4 


x (j--m—2)}L( +2), 
ef gee ie! — ni | 


15 
xv {( Ry cman OSS 


ens . fe ae 
(j,m|E|j—1, m) =——FE —_ =n | /{(j—m)(j+m)}L (0). 


266 On the Matrix Elements for Octupole Radiative Transitions 
Aj=—2 
(j,m|E|j—2, m43)=}F vy (jAmt+ Gm) jFm—1)\(jFm—2) 
x (JF m—3)(jFm—4)}L( +3), 
(j,m|E|j—2, m+2)= 44K (2j£3m+2)V/{(jFm)(jFm—1)(jFm—2) 
x (JF m—3)}L(+2), 
(j,m|E|j—2, m£1)=4fF(j43m+1)V/{(j+m)(j—m)(GFm—1) 
x (JF m—2)}L( +1), 
(j, m|E|j—2,m) =—3Fmy {(j+m)(j—m)(j-+-m—1)(j—m—1)}L (0). 


Aj=—8 
(j,m|E|j—3, m+3)= +4G/{(jFm)(j—-m—1)(jFm —2)(j-m—3) 
x (j-Fm—4)(jFm—5)}L( +3), 
(j, m|E|j—3, m+2)=}G/{(j+-m)(j —m)(jFm—1)(j=m—2)(jFm—3) 
x (j--m—4)}L(+2), 
(j,m|E|j-3, mt 1)=+3GV/{(j+m)(j—m)(j+m —1)(j—m 1) 
X (9 m—2)(j-m—3)}L( +1), 
(j, m|E|j—3,m) =—Gy/{(j+m)(j—m)(j-+-m—1)(j—m -1) 
xX (j--m —2)(7—m—2)}L(0). 
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XXX. Fluctuations in the I ntensity of Radio Waves from Galactic Sources. 


By C. G. Lirrie and A. MAXWELL, 
Jodrell Bank Experimental Station, University of Manchester *. 


[Received January 18, 1951.] 


SUMMARY. 


It has been shown previously that the fluctuations in the intensity of 
radio energy from the galactic sources have a local terrestrial origin. 
This paper describes further observations of the two most intense sources, 
in the constellations of Cygnus and Cassiopeia, taken over a wide range 
of angles of elevation. At high angles of elevation the radio fluctuations 
are shown to be correlated with the occurrence of “ spread ’”’ ionospheric 
echoes from the F region. When the sources are low on the northern 
horizon fluctuations are always observed; these are probably introduced 
by the passage of radio waves through the continuously disturbed 
ionospheric regions at high magnetic latitudes. Spaced receiver obser- 
vations taken over base lines of 0:1, 4 and 11 km. enable the scale of the 
radio energy diffraction pattern across the ground to be determined. 


§ 1. InTRODUCTION. 


’ FLUCTUATIONS in the galactic radio emissions from the region of Cygnus 
were first observed by Hey, Parsons and Phillips (1946). They concluded 
that these fluctuations originated in a high intensity variable source of 
small angular diameter. Subsequently Bolton and Stanley (1948), 
using an interference technique, discovered a localized source of angular 
width <8’ in the Cygnus region, from which they recorded a variable 
component of radiation at 100 Mc./s. At this time it was generally 
assumed that the variable component was due to true variations in the 
‘source emission. Recent combined work by Smith (1950), Little and 
Lovell (1950), has shown, however, that the fluctuations in the radio 
emission from a galactic source are not inherent in the source itself, but 
are impressed on the incoming radiation as it passes through the earth's 
ionosphere. 

The occurrence of fluctuations at night when the electron density of the 
lower ionospheric layers is small suggests that the F region is the most 
probable origin of the fluctuations. (The existence of fluctuations when a 
source is high in elevation indicates that tropospheric refraction is not the 
cause of the phenomenon.) It was suggested by Little and Lovell that 
diffraction of the incident radiation by irregularities in the F region 


——— 


* Communicated by Dr. A. C. B, Lovell, 
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would suffice to account for the observed fluctuations. An analysis of 
the occurrence of fluctuations and of “spread ”’ F echoes which are 
believed to be associated with such irregularities has now been completed 
and is given in § 4. 

The spaced receiver observations referred to above were made at transit of 
the two most intense galactic sources, in Cygnus (Declination approximately 
+40° 30’, Right Ascension 19" 57m) and in Cassiopeia (Declination 
approximately +58°, Right Ascension 23 21™), This work has now . 
been extended to all angles of elevation by using aerials movable in azimuth 
and elevation. In this way it has been possible to investigate in more 
detail the effect of the ionosphere on the radiation. 


§ 2. APPARATUS. 


The apparatus consisted of two equipments working on a frequency 
of 81-5 Mc./s., each using a directional aerial, movable in azimuth and 
elevation, a sensitive, high gain receiver and a recording meter. One 
receiver was sited permanently at the Jodrell Bank Experimental Station, 
Cheshire (Lat. 53° 14’ N., Long. 2° 18’ W.), and was connected to a 30 ft. 
aperture, focal plane paraboloid aerial (free space beam width +11°, 
+15° to half power). A mobile equipment with an array of two Yagi 
aerials (beam width +11° +21°) was used at spacings up to 11 km. 
from Jodrell Bank. 

The receivers used were of the conventional superheterodyne type, 
modified to include two preamplification stages of low noise factor, and 
two stages of post detection amplifcation for recording weak signals. 
The characteristics of the receivers were: Frequency 81-5 Me./s., band- 
width 0-6 Me./s., output time constant 0-2 seconds, noise factor 4, overall 
gain 10°. All receiver power supplies were carefully stabilized. Pen 
recording milliammeters were used as recording instruments, with chart 
speeds ranging from 12 inches per minute to 3 inches per hour. 

In this paper, when it is stated that galactic noise fluctuations were 
“not observed’ it is implied that the fluctuations were less than the 
ripple in the receiver noise output. For the apparatus used, the minimum 
detectable level was of the order of 1/500 of the total receiver noise, 
corresponding to a 5-10 per cent modulation of the source intensity. 
Strong galactic noise fluctuations showed peak intensities up to 40 times 
the minimum detectable fluctuation level. 


§ 3. MetTHop OF OBSERVATION. 


The experiments were carried out over the period May 1949—December 
1950, and may be divided into 3 sections: (i.) Transit observations to 
investigate the correlation between the fluctuations and local ionospheric 
phenomena, (ii.) Continuous observations to investigate the changes in 
amplitude and period of the fluctuations with elevation of the source. 
(iii.) Spaced receiver observations to determine the linear dimensions of 
the ionospheric irregularities responsible for the fluctuations, | 
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§ 4. OBSERVATIONS aT TRANSIT —CoRRELATION WITH SPRBAD F. 


The experiments of Smith, Little and Lovell, referred to in $1, showed 
that the fluctuations must originate in a local medium. An analysis of a 
series of transit observations, taken over the period May 1949—December 
1950, has now been made to establish the extent of the correlation of the 
fluctuations with the ionospheric discontinuities which give rise to a 
“ spread ” in the critical frequency and height of F region echoes. 

Table I, compares the occurrence of galactic noise fluctuations and 
spread F conditions during the transit of the Cygnus and Cassiopeia 
sources. The galactic noise observations, made at Jodrell Bank, refer to 
the period +1 hour from transit, whilst the ionospheric observations 
were made at Slough, some 230 km. to the South, and refer to the period 
transit +2 hours. 


TABLE I, 


Correlation of galactic noise fluctuations with spread F. 


No. of No. of occasions Spread F Galactic 
occasions when| when spread F and observed noise 
Radio spread F and noise fluctuations but no fluctuations 
source | galactic noise were galactic observed but 
observations | (i.) Both | (ii.) Both noise no 


were made | observed| absent | fluctuations | Spread F 


Cygnus 205 53 (26%) | 106 (52%)! 16 (8%) 30 (14%) 


Cassiopeia 195 69 (35%) | 74 (38%) 36 (19%) 16 (8%) 
Transit observations May 1949 to December 1950. 


It is seen that the galactic noise fluctuations and spread F were 
simultaneously present, or absent, for 75 per cent of the observations, 
i. e. only 25 per cent of the runs indicate the presence of the one without 
the other. (Smith, Little and Lovell have shown that there is a 10 per 
cent chance that fluctuations will be observed at only one of two sites 
210.km. apart and hence a correlation figure higher than 90 per cent 
cannot be expected.) ; 


§5. Spaced RECEIVER EXPERIMENTS. 


Previous work with spaced receivers showed that there was no correlation 
between (transit) fluctuations observed on two receivers separated by a 
distance of 210 km. but that a fairly good correlation existed for a spacing 
of 4km. These experiments have been continued and extended to all 
angles of elevation to discover where correlation ceases. The results may 
be summarized as follows : 

(i.) Two receivers located at the same site, observing the same source, 
record identical fluctuations. 
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(ii.) For a separation of 110m., the correlation coefficient is still 
extremely high (>-98). The deviation from unity is within the range of 
the experimental limitations. 

(iii.) As already reported (Little and Lovell 1950), at 4 km. the correla- 
tion between the fluctuations is not complete, but remains high with a 
correlation factor varying between 0-5 and 0-95. 

(iv.) Over an 11-2km. North-South base line there is no correlation 
between transit fluctuations. However, when a source is near lower 
culmination the “effective” base line (i.e. the projection of the 
actual base line on a plane normal to the line of sight) becomes small, 


Fig. 1. 


ol 25™ olh 30” UE Classe ol"40™ 
1950 March 16 


rego) NCO} 20° 45™ U7 20 Som 20°55" 
1950 March I5 
Cassiopeia fluctuation records taken simultaneously by two receivers 11 km. 
apart on a North-South base line. 
Records A: Correlation 0-9 (effective base line 4-1 km.). 
Records B : No correlation (effective base line 7-7 km.). 


and the correlation factor may rise to as much as 0:95. This is illustrated 
in fig. 1, in which correlations between fluctuations for different effective 
base lines are compared, and in Table II., which shows the change in the 
correlation as the Cassiopeia source nears lower culmination. (The base 
line used was actually 10° KE. of North so that the effective base line is 
given by 11-2 sin [cos~! (cos « cos B—10)], where « is the angle of elevation 
and £ the azimuth angle.) ) 

(v.) At 210km., between Cambridge and Jodrell Bank, there is no 
correlation between transit fluctuations, 
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TABLE II. 


Change of cross-correlation coefficient between fluctuation records 
with change of effective base line. 


Effective 
base line 
Kilometres 


Elevation | Azimuth 
degrees | degrees 


Correlation between 
fluctuation records 


11-2-6-4 correlation not significant 
6-4 0- 


oo 


5 
4-4 
4-] 
4-6 
5: 
6- 


correlation not significant 
correlation not significant 


Hoge 


§ 6. CONTINUOUS OBSERVATIONS. 


During the year November 1949-November 1950 some 3000 hours of 
observations were made on the Cygnus and Cassiopeia point sources by 
following them with steerable aerials. These observations usually took 
the form of a 24-hour run on each source once or twice per week throughout 
the year. Daytime observations were limited during some months by 
strong radio emissions from active areas on the Sun (of the order of 103 
intensity of the galactic source radiation), by atmospherics from thunder- 
storms during the summer months, and occasionally by strong -local 
interference. 

The results of this series of observations may be summarized as 
follows :— ; 

(i.) For a source at high angles of elevation the fluctuations are very 
closely related with the presence of spread F, and as this latter is essentially 
a night time phenomenon, fluctuations at high angles of elevation are 
usually confined to the hours of darkness. _ Thus, near upper culmination 
fluctuations were observed on 35 per cent occasions when transit occurred 
during the hours 1800-0600, but on less than 5 per cent for transit during 
i hours 0600-1800. 

4 eee angles of elevation when the line of sight crosses the 
ionosphere at high magnetic latitudes, fluctuations were always observed, 
independent of the hour of day. These fluctuations menos general less 
marked in the case of Cassiopeia (minimum elevation 21°) than for Cygnus 
ini ion 4°). 
Leanne Split Fe the fluctuations of the C ygnus source keprractod 
- for aerial ground lobes) plotted against the elevation e fine perce is 
shown in fig. 2 (a). It is seen that for angles below 20 the average 
amplitude of the fluctuations, expressed as a fraction of the mean Lee y 
of the radiation from the source, increases rapidly, reaching a maximum 
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of about 90 per cent at lower culmination. It will also be seen that for 
lower culmination in daylight the maximum amplitude of the fluctuations 
is significantly less than for lower culmination at night. ; 

(iv.) The marked elevation effect associated with fluctuations from the 
Cygnus source is further illustrated in fig. 3. In this diagram the “ onset ~ 
and “ cessation ”’ of the fluctuations (taken here as the time at which the 
fluctuations respectively increase or decrease to an amplitude of 25 per 
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(Upper diagram). Average amplitude of Cygnus fluctuations plotted against 
elevation of the source. 


(Lower diagram), Average number of Cygnus fluctuations per minute plotted 
against elevation of the source. 

In each case curve A represents the average values during the six winter 
months (October—March), and curve B the average values during the 
six summer months (April-September). Cygnus lower culmination occurs 
in daylight during the summer, and in darkness during the winter. 


cent) is plotted throughout the year. It is seen that the onset of the 
fluctuations occurs at approximately the same sidereal time each day, 
about 5 hours before lower culmination, corresponding to an elevation 
of 23°, and that the fluctuations cease approximately 5 hours after lower 
culmination. The spread of points during the August-November period 
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on the onset curve, and during March-June on the cessation curve, 
correspond to fluctuations at high angles of elevation during the hours of 
darkness. 

(v.) The duration of individual fluctuations at any given angle of 
elevation and azimuth varies considerably from night to night. For 
example, individual fluctuations at transit of the sources may last for 
2 or 3 minutes on some nights, whilst on other nights for only 
20 seconds. Typical durations, however, range from about 1 minute at 
high angles of elevation to about 15 seconds at low angles of elevation. 
The number of fluctuations of the Cygnus source per minute plotted as a 
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‘Times of ‘‘ onset’ and “ cessation’ of Cygnus fluctuations. The two dotted 
lines correspond to a source elevation of 23° (7. e. 5 hours before and after 
lower culmination respectively). 


function of the angle of elevation is shown in fig. 2 (b) ; it will be seen 
that the number of fluctuations per minute for lower culmination in 
daylight is somewhat less than for lower culmination at night. 


§ 7. DISCUSSION. 
(i.) The Mechanism of Fluctuations : Diffraction im the Ny onosphere. 
It was suggested by Little and Lovell that the galactic noise fluctuations 
could be caused by diffraction of the incident radiation by localized 
irregularities in the refractive index of an ionospheric layer. ete 
irregularities would cause a distortion of the incident wavefront due to t 1€ 
different optical path lengths in the ionosphere, and the resultant cisirsedon 
pattern would be observed as an irregular distribution of amplitude across 
the ground. This distribution would not be constant with time,—due to 
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the presence of winds and turbulence in the ionosphere, and to the apparent 
change in position of the star as the Earth rotates,—and hence a fixed 
receiver would observe fluctuations in the intensity of the incident 
radiation. 

The detailed implications of these ideas have been discussed by Little 
(1951) on the basis of a theorem by Booker, Ratcliffe and Shinn (1950), in 
which it is shown that the auto-correlation function of the diffraction 
pattern across the ground is the same as that of the distortions of the wave- 
front emerging from the ionosphere ; that is, the scale of the pattern at 
ground level is the same as that of the distortions in the emergent wavefront. 
On the assumption that the scale of the distortions in the wavefront is the 
same as that of the irregularities in the 3-dimensional ionospheric screen, 
then the experimental results described in this paper will give some 
indication of the lateral dimensions of the irregularities. 


(ii.) Correlation between records from spaced receivers. 

In §5 it is shown that the correlation between fluctuation records begins 
to fall off at spacings of approximately 4km., and it would appear, 
therefore, that the ionospheric irregularities causing the fluctuations have 
lateral dimensions of the order of 4 km. 

If we now compare the observations made at the transit and the lower 
culmination of, say, Cassiopeia by two receivers 11-2 km. apart on a 
North-South base line, then allowance must be made for the changes in the 
* effective’’ base line (§5 (iv.)). Thus at transit, when the source is at an 
elevation of 85° and the effective base line is nearly 11-2 km., we should 
expect poor correlation if the F region irregularities are of diameter 
approximately 4km. Towards lower culmination, however, when the 
source is at an elevation of 21°, and the effective base line tends toward a 
minimum of 4-1 km., the correlation would be expected to be good. This 
is in agreement with the experimental observations. 

The above correlation results indicate that the scale of the diffraction 
pattern, in a plane normal to the incident radiation, is independent of the 
elevation of the source for elevations greater than 20°. In the case of 
Cygnus lower culmination (elevation 4°), the sight line crosses the iono- 
sphere at high magnetic latitudes, and the correlation only becomes 
significant when the effective base line is | km. or less. This suggests that 
the ionospheric irregularities at these high magnetic latitudes are appre- 
ciably smaller than those at medium latitudes. 

(iii.) Fluctuations at low angles of elevation. 

There are clearly two important differences between transit observations 
and those taken near lower culmination : (a) as shown in fig. 4, the sight line 
from Jodrell Bank to the Cassiopeia and Cygnus sources at lower culmination 
crosses the F region of the ionosphere at latitudes of 62° N.and 71° N. respec- 
tively,—the latter well within the auroral zone. At these high latitudes the 
ionospheric structure is not well defined and fluctuates rapidly with time ; 
sporadic F is much more intense and occurs during the day as well as:at 
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night ; sporadic E densities are also greatly in excess of the densities at 
lower magnetic latitudes.. (6) When the sources are low in elevation the 
depth of ionosphere along the line of sight is increased. The variation of 
the amplitude of the Cygnus fluctuations with the change in the effective 
thickness of the ionosphere is shown in fig. 5. (This latter jeauneen 
plotted for a disturbing region of two assumed heights and thicknesses : 
viz., between 300 and 500 km. and between 250 and 1000 km. : owing 
to the curvature of the earth the thickness factor does aN follow 
a simple cosecant law.) It is seen that the effective thickness and 
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Diagram showing positions at which radiation from the Cygnus and Cassiopeia 
sources incident at Jodrell Bank crosses the F region of the ionosphere 
(assumed at a height of 400km.). At low angles of elevation the sight 
lines cross the F region near the auroral zone. 


amplitude curves follow each other closely for elevations greater than 20°, 
and it would therefore appear that down to an amplitude of about 20° the 
increase in fluctuation amplitude could be accounted for in terms of 
increased thickness of the disturbing region. Below 20°, however, where 
the rate of increase in the amplitude of the fluctuations is much greater 
than the corresponding rate of increase in effective thickness, it would 
seem that the disturbed nature of the ionosphere at high magnetic 
latitudes is the main contributing factor. 


276 C. G. Little and A. Maxwell on the Fluctuations 


The effect of the troposphere at these low angles of elevation is not 
readily determinable. However, on metre wavelengths, and for a 
minimum elevation of 4°, it is unlikely that fluctuations of the type 
observed could be introduced by the troposphere, because (a) the variations 
in the refractive index of the troposphere are too small: of the order of 
10-® (cf. Booker and Gordon 1950), compared with 10~4-10~ for the 
ionosphere, and () the effective thickness of the troposphere traversed at 
4° elevation is about 100 km., compared with 1000 km. for the effective 
thickness of the F region of the ionosphere. 


Fig. 5. 
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Curve A: Amplitude of Cygnus fluctuations with elevation. source setting, 
® source rising (averages over 12 months). 

Curves B and C: Ratio of effective thickness of disturbing region in line of 
sight to zenith thickness—(Curve B for a disturbing region of assumed 
minimum height 300 km. and maximum height 500 km.; curve © for 
assumed minimum height 250 km. and maximum height 1000 km.). 


(iv.) The duration of individual fluctuations. 

Time variations in the intensity of the galactic radiation could be 
brought about in several ways: (a) a random growth and decay of the 
irregularities, (b) a turbulent motion of the irregularities, (c) the movement 
of the ionospheric irregularities as a whole under the influence of an 
ionospheric wind, (d) the apparent motion of the star relative to the 
ionospheric irregularities due to the rotation of the earth. 


—_ ——— 
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Tonospheric observations have shown that changes in the electron density 
of the F region are of periods very much longer than 30 seconds, the average 
duration of an individual fluctuation. This, together with the observation 
that correlation between fluctuation records is high over base lines of 
4km., suggests that fluctuations are not introduced to any great extent 
by turbulence, or by random growth and decay. 

The duration of individual fluctuations should therefore be determined 
by the ratio of the size of phase irregularities in the ionosphere to the 
components of wind and star velocity in their plane. Wind velocities 
of approximately 7 km./min. have been found by Munro (1948), and for 
4 km. irregularities would lead to fluctuations with a period of the order 
of 30 seconds. In the absence of winds, the apparent motion of the star 
due to the rotation of the earth would lead to fluctuations with a period of 
approximately 3 minutes. These figures agree fairly well with the 
observed limits of from about 20 seconds to 5 minutes for the fluctuations 
at transit. 

The results described in § 6 (v.) and illustrated in fig. 2 (b) show that 
the period of the fluctuations is significantly less near lower culmination 
(decreasing steadily for elevations less than about 20°) than at transit. 
This effect is to be expected for three reasons: (a) the radiation then 
traverses the ionosphere near the disturbed auroral regions, (b) the depth 
of the ionosphere along the line of sight is considerably increased, and 
(c) ionospheric winds of 20 km./min. are reported to exist at high latitudes 
(Meek 1949). 

§ 8. CONCLUSION. 

The main conclusion to be drawn from the experiments described above 
is that the fluctuations in intensity of the radio emission from the galactic 
sources are introduced when the radio waves passthrough an inhomogeneous 
ionospheric region, and spaced receiver experiments suggest that these 
inhomogeneities are of the order of 4km. in size. Correlation with 
ionospheric phenomena indicates that the F region is mainly responsible, 
but it is not possible to exclude other ionospheric effects in cases where 
the radiation passes through disturbed auroral regions. 

~ Further investigation of galactic noise fluctuations on various frequencies, 
and over varying base lines, should yield valuable information on the 
structure of the upper regions of the ionosphere. 
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XXXI. The Electric Quadrupole Moment and the Magnetic Moment of 
the Bromine Nucleus. 


By J. D. Ranape*, 
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[Received January 9, 1951.] 


ABSTRACT. 


Deviations from the interval rule have been observed in the hyperfine 
structure of the bromine are spectrum. These are attributable to an 
electrical quadrupole moment of the bromine nucleus. Coupling constants 
are determined for the term involved. The nuclear quadrupole moment 
is calculated to be 0-28 x 10-*4 cm.2, and the magnetic moment 1-8 nuclear 
magnetons. 


$1. INTRODUCTION. 

DEviATIONS from the interval rule in the BrI spectrum were first observed 
in the term 5s*P;,. by Tolansky and Trivedi (1940). These can be 
attributed to the nuclear electrical quadrupole moment. The separations 
in this term were measured accurately from lines in the red and the 
infra-red regions where the resolving power of the Fabry—Perot interfero- 
meter is high. There is no evidence of any possible perturbation from 
another configuration. Table I. shows the observed positions of the 
hyperfine structure levels as compared with those calculated from the 
interval rule. 


TABLE I. 
ev alllcanecte:. 4 33 2 i 
Observed...... 417 218 83 0 
Caleculated..... 417 232, 93 0 


(Interval rule) 


At the time the quadrupole moment was not evaluated for want of 
the coupling constants of the term. Hence, here, following a method 
similar to that used by Schmidt (1939) for iodine, coupling constants 
of the term involved and thus the quadrupole moment are determined. 


§2. QUADRUPOLE MOMENT. 
The interaction between the angular momentum J of a non-spherically 
symmetrical configuration and the spin I of a non-spherically symmetrical 
nucleus takes the form 


E=@) + He | be(c+1), 


where 


c=F(F+1)—1(1-+1)—J(S-+1). 


* At present, at the Institute of Science, Bombay, India. 
+ Communicated by Professor 5. Tolansky. 
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In this interaction law ay, a, and b are constants, a) being the displacement 
of the centre of gravity of the hyperfine structure pattern from the 
position the term would have occupied if the spin had been zero, a, is 
the hyperfine structure interval factor which is a function of the nature 
of the electron configuration and is also proportional to the nuclear 
magnetic moment, b is a measure of the deviation from the interval rule. 
The quantity b is a function of the nuclear electrical quadrupole moment 
of the nucleus. 

The values of the interaction constants given by the observed intervals 
are (units cm.—1 x 10-3) 


ay= 223, 
a,= 47-1, 


For the calculation of the quadrupole and the magnetic moments, it 
is necessary to know the coupling constants for the configuration electrons. 
For these, first, the coupling constants of the 4p+ electron group in the 
singly ionized bromine atom is calculated, then the attachment of a 
5s electron to this group, giving rise to the 4p*. 5s configuration of the 
neutral atom is considered. It is assumed that the attachment of a 
5s electron to the 4p* group does not alter the coupling constants in the 
group. The p* electron group can be treated as equivalent to the p? group. 
Hence for calculation p*.s can be treated as p?.s. The effect of this on 
the quadrupole moment formula is a change of sign since the spin-orbit 
interaction changes sign. The eigenfunctions for the p*s configuration 
are written in terms of the single electron functions. The term 4P;). 
arises from the term ?P, of BrII. Thus for J=2 and m=2 the eigen- 
functions for the configuration p? are, as given by Breit and Wills (1933) 


pe 
ties S/22"1 [20 i 1/9 eae 
E Oe Pseorr (P33P ‘1/2 —Pi/eP 3/2), 


33]? 1 3/2 1/2__ 1/2132 
55) | ae V3 (p3/2P 3/2——P3/2P 3/2). 


9 
a 


Here p3/5 etc. signify the single electron eigenfunctions, the magnetic 
quantum number being given by the superscript to the right. The 


eigenfunction of the °P, term is a linear combination of the above functions 


and is given by 
5 fe 33]? 
3p2) — wees Be pes 
vewd-o[aa] vLaa]. 
The coefficients c, and c, are to be determined from the multiplet structure 
and from the value of the spin-orbit interaction. The values of the 
spin-orbit interaction (6) and the electrostatic energy integral (F,) are 
given by Robinson and Shortley (1937). Using the values 
6=2958cm.-? and F,=1750 cm.-} 
and the term values for BrII, leads to 


¢;=0:434 and c,=0-901, 
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The eigenfunctions of the *P5/2 term arising from the attachment of the 
s-electron to the p? group are given by 


HAPs) —=si/2 - #PP3). 
The formula for the electrical quadrupole moment q as given by Casimir 
(1935) for a p-electron is 


. b.Z,.H 2i(2i—1)j(2) 1) x 0-98 ‘ 
ee ee = i = ois =e 

§(3 cos? d—1) Ie MRIS ald 
where H is a relativity correction. The average value of (3 cos? d—1) 
is determined by using the above eigenfunction in conjunction with 
formule given by Schiiler and Schmidt (1936). For 4P,, we have 


(3 cos” d—1)=—2(2/2c,c,8,—R;), 

where S, and R, are the relativity corrections and c, and c, are the 
constants already computed. Thus finally 

__ bZ,H2i(2i—1)j(2j—1) x 0-986 x 10-4 

8, —3(29/2¢,¢,8, —c{R:)] 
With c,=0-434, c,=0-901, S=1-07, Rj=1-045, H=1-024, i=3, j=3, 
Z;=31 [(Z—4) for a p-electron] and 6=2958 cm.—! giving 
G07 1612) 10-24.cm.? 


and with b=0-17 
q=+0-28 x 10-74 em.?, 


No separation of the isotopes is revealed by hyperfine structure. The 
value obtained here compares well with the values found by the micro-wave 
technique and reported by Townes, Holden and Merrit (1948), as 


Brae g—-.0 25% 10-* em.*, 


Brey. G= 4-023 x 10-74 cm”. 


§3. Magnetic MoMENT. 

‘Using the earlier data of Tolansky (1932) for the separations of the 
terms !P5)9, *Ps/2 and ?P3). of Bri, Schmidt (1938) calculated the nuclear 
magnetic moment of bromine as 2-6 nuclear magnetons. Since then 
Tolansky and Trivedi (1940) have remeasured and redetermined the 
separations more accurately. The interval factor for the term *P,), 
is now given to be 20 cm.-1x 10-3 instead of the value 42 ern! Sal Os 
first given and used by Schmidt (1936). Therefore the nuclear magnetic 
moment has been re-determined using the improved data now available. 
The total splittings for the terms *P;,, and *P;,, which are now being used 


are 
(Py 4)=417 cm, X10, 
and 
(2P 273) ==120 cm, X 108% 
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bo 


Breit and Wills (1933) have calculated the interval coupling factors for 
the #P,). and 4P;,. terms, in intermediate coupling, giving 


A(4P 5/o)=ta(s)-+2[(3c?-+403)a’ +-c2a" +4-V/2c,c90'""], 
A(4P5/2)= —3a(s s)+2[ (8) + 4e3)a a" +a” +4 2c,0,0'"]. 


Here a(s) is the interval constant for the s-electron and a’, a’’ and a’”’ are 
the interval constants for the p-electrons, c, and c, are the coupling 
constants as computed above. The constants a’, a’’ and a’” are related 
as follows. 
a’ =a’ ; 5F’/E’ 

and a’ =—a' . 5/16. G/F’. 
Here F’, F” and G are the relativity corrections. 

Using the values c,=0-434, c,=0-901, F’”=1-09 and F’=G=1-02 leads 
to the interval factors 


A(4P.,)9) =0-2a(s)-+0-823a’, 
A(4P3)9)= —0-2a(s)+-1:235a’. 
From these relations and the known values of 
A(4P./2.)=47-1 om.-?X10-% and A(4P3).)=20 em.-? x 10-3 
given by Tolansky and Trivedi, it follows that 
als)=101-3 ems 4 1027; 
@ ==32-6 om. 1x 10: 
For a(s), Goudsmit’s formula leads to 
a(s).0. 0.6. 0-1185 
V= 77381 (1—deldn)a, we 
For a’, 


a’ .i. ZH. 3-446 
oF eda Om [a a . . . . . . . (2) 


In equation (1), % is the effective quantum number for the s-electron 
when it is the only valence electron present. This happens in BrV. 
The corresponding configuration is 4s*.5s. In neutral bromine the 
configuration is 4s. 4p4.5s. The absolute value of this term of BrV 
is not experimentally known. It is, however, here determined by 
extrapolation using the known term values for the 5s? . 5s, 5s?.. 6s ete. 
terms of the iso-electronic sequence gallium I, germanium II and 
arsenic IIT. The term values for these are taken from Goudsmit and 
Bacher (1932). From this ,_~=3:95 and (1—ds/dn)=-0-65. With i—3, 
Z=35, Z;=31, H=1-024, F’=1-02, F’=1-09, m4q=3-95, (1 —ds/dn)=0-65, 
and Z)=5 (Br V), giving 
from a(s): =1:79, 


and from a’: ~=1-77, 
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“a 


The agreement between these values is perhaps fortuitous and much 
better than can be expected. One can take closely enough 


/4=1-8 nuclear magnetons. 


The determination of the nuclear magnetic moment of bromine has been 
carried out by the nuclear paramagnetic resonance absorption method 
and the magnetic beam resonance in CsBr and LiBr. The values given 
by Pound (1947) are 


(Brs')=—=2-2 nuclear magnetons, 
(Br?®)=2-1 nuclear magnetons. 


When the approximations used here are taken into account the difference 
is not significant and indeed the two methods give reasonably concordant, 
results. 
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ABSTRACT. 


Measurements for the hyperfine structures of 18 classified and 21 
unclassified lines of the first spark spectrum of bromine are reported. 
The spectrum was excited by a high-frequency electrodeless discharge 
in pure bromine vapour. The patterns are analysed and interval factors 
for fourteen terms are given. 


§1. INTRODUCTION. 

MrasureMeEnts for the hyperfine structures of several lines of the arc 
spectrum of bromine have been reported earlier by Tolansky (1932) and 
Tolansky and Trivedi (1940), who proved that both the bromine isotopes 
79 and 81 have the same nuclear spin, 3/2, and have practically identical 
nuclear magnetic moments. Tolansky and Trivedi also reported measure- 
ments of the hyperfine structure for four lines of the first spark spectrum 
of bromine. 

The measurements reported here consist of a study of the classified 
and the unclassified lines of the spark spectrum of bromine. The first 
spark spectrum of bromine was partially classified by Lacroute (1935), 
who studied the Zeeman effect but could not sare, them fully due to 
lack of hyperfine structure data. 


$2. EXPERIMENTAL. 


The spark spectrum was excited in bromine vapour with an electrodeless 
high-frequency discharge which is gentle and can be considered essentially 
as a low voltage excitation in which arc lines tend to predominate. A 
brilliant arc spectrum can be excited in bromine vapour at a pressure of 
about | mm. As the gas pressure is lowered, the increase in the mean 
free path enables the ions to acquire higher energies between collisions 
and spark lines appear. Sharper spark lines result than those excited 
in a Geissler discharge tube. By reducing the bromine vapour pressure 
below | mm. and by using a capillary more constricted than that necessary 
for an are spectrum, strong spark lines were excited. 

Oscillations were generated by a 500 watt thermionic valve, the 
wavelength being of the order of 20 m. The field was applied by means 
of an external electrode of copper foil to a pyrex discharge tube which 


*Communicated by Professor. 8. Tolansky. 
} At present at the Institute of Science, Bombay, India, 
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was incorporated in a vacuum circulating system having charcoal and 
liquid oxygen traps in the circuit. A continuous stream of pure bromine 
vapour was circulated through the discharge tube from a side limb 
cooled by immersion in a flask containing liquid oxygen. By adjusting the 
height of the bottom of the limb above the liquid oxygen surface, the bromine 
could be cooled to a temperature of —80° C. at which a good excitation. 
of the spark spectrum resulted. Cooling to —45° C. is necessary for the 
are spectrum. No carrier gas was used. The portion of the discharge 
tube emitting the spectrum was cooled by directing on to it a strong 
blast of air. 

The high resolving power instruments used were a variable gap Fabry- 
Perot interferometer coated either with silver, aluminium or an 
aluminium-magnesium alloy, and a quartz Lummer plate by Hilger, 
(20 cm. long and 3.42 mm. thick). The interferometer was used in an 
external beam mounting crossed with a large Hilger glass-quartz Littrow 
spectrograph. 

§3. OBSERVATIONS. 

Table I. contains hyperfine structures measured for classified lines of 
BrII spectrum. The wavelengths and allocations are taken from Lacroute. 
The units used are cm.~! x 10-3 and beneath each component is indicated 
a visual estimate of its intensity. The lines are classed A, B, C according 
to whether the resolution in the pattern is good (A), moderate (B) or 
poor (C). The factors affecting the class are the complexity and the 
narrowness of the patterns and the widths of the lines. For lines classed 
as A, the errors are some | or 2 cm.~1X10-%. Whenever two components 
are not completely resolved, the blend is indicated by binding of 
components with a hyphen. Measurements for a number of unclassified 
lines are given in Table II. 


$4. ANALYSIS. 


For the singly ionized bromine atom Hund’s theory predicts three 
families of terms arising from the basic terms 48, 2D and 2P of BrlII. 
All the lines classified by Lacroute (1935) belong to the 48 system. The 
analysis of lines involving common terms will now be considered. 

Of the three lines A 5164-4, 5238-2 and 5182-4 one has a term with a zero J 
value and hence no structure. The three lines have a common lower 
term 5s38,, the upper terms being 5p°Po, °P;, *P, respectively. The 
structure of the line \ 5164-4 5s3S,—5p?P, gives that of the term 5s°8,. 
From this, the structures of the upper 5p°P, and 5p*P, terms follow as 
shown by fig. 1. The agreement between the predicted and the observed 
patterns is reasonably good and 5.8, is not fully resolved. 

The term 5358, will give rise to regular quartets in lines for which the 
upper terms have narrow structure. The two lines A 4193-5 and 4230-0 
have a common lower term 535S,. Since the structures of the upper 
terms 5pP, and 57°P, have been derived above, these lines give the 
structure of the term 55°S,. The two central components of A 4230 are 
not ‘properly resolved. : 


SER. 7, VOL. 42, NO. 320.—MAROH 1951 x 


286 


J. D. Ranade on the 


The lines \ 4816-7, 4785-5 and 4704-9 have well resolved quartet patterns. 
The structures in these lines are essentially those of the term 5s°S, and 
The analysis is 


the number of components directly gives the spin. 


TABLE I. 


Hyperfine structures, of classified Br II lines. 
Unitscnm it. 


Wavelength Allocation. 
5238-2 5s°8,—5p®P, 
5182-4 5s°8,—5p*®P, 
5164-4 5898, —5p>P, 
4848-8 5p®P,—6s°S, 
4816-7 58°S.—5p*P, 
4785-5 5s°S,—5p°P, 
4766-0 5s°P,—6s°S, 
4735-4 5p®P, —6s5S, 
4704-9 5s°8,—5p*P, 
4396-4 5p®P,—4d®D, 
4356-9 5p>P,—4d°D, 
4230-0 58°S,—5p®P, 
4193-7 5s°S,—5p'P, 
3980-4 5p®P,—4d°D, 
3980-0 5p°P,—4d°D, 
3935-2 5p®P,—4D, 
3914-5 5p'P,—40D, 
3914-1 5p®P, —4d5D, 


Structure 
Red 
0 95 
(1) (2) 
0 38 93 
(1) (2) (3) 
0 48 
(1) (2) 
Borad single-width 
50 
0 154 264 
(5) (3) (2) 
Q).. “125° 9s213. 
(4) (3) (2) 
0 65 
(2) (1) 
0 76 
(2) (3) 
OPP 130 e218 
(5) (3) (2) 
0— 75 
(1) (1) 
Single 
0 188 270 
(Sit (2) at 
0 130, -219 
(4) (3) (2) 
0— 49 
(2) (1) 
Single 


Width 50 
Single 


Single 


Violet 


Class 


35 eel fen Aigsh Maes 


Q 


illustrated in fig. 2 giving the structures of the upper terms 5p5P,. 5P, and 
°*P;. The agreement between the predicted and the observed values is 


good, 
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TABLE IT. 


Hyperfine structures of unclassified BrII lines. 
Units cm‘=? x 10-8, 
ae ee Se er et 1 


Structure 
Wavelength Red Violet 
4921-3 0 112 
(1) (3) 
4802:3 434. 337, 188 0 
(1) (2) (3) (5) 
4742-7 0 122 198 237 
(5) (3) (1) (1) 
4728-2 Single 
4719-8 0 59 119 
(2) (2) (5) 
4693:3 0 144 218 384 
(6) (3) (1) (1) 
4678-7 Single 
4542-9 0 121 254 362 
(1) (2) (5) (9) 
4365-6 0 67 156 
(2) (2) (3) 
4351-2 0 176 344 
(5) (3) (1) 
4297-8 0 156 246 
(3) (1) (1) 
4291-4 0 158 266 364 
(9) (6) (4) (3) 
4223-9 118 0 215 341 
(3) (10) (5) (2) 
4179-6 0 146 232 
(7) (2) (5) 
4140-4 Single 
4135-9 Single 
3986-3 132 0 205 338 
(1) (10) (3) (5) 
3970-6 Single 
3968-3 0 129 224 
(9) (6) (4) 
3955-4 Single 
3950-6 Single 


The lines \ 4766-0, 4735-4 and 4848-8 have complex structures. Except 
for \ 4735-4, the components are closely packed so that resolution is poor. 
4848, appears as a broad single line pattern of width 50 em.-! x 10-%, 
The structure of the 6s5S, term is obtained from these lines, 


xX 2 
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The group of 445D terms combines with the 5p* Pore and 5p°Pyio3 
multiplets which have narrow structures. Most of the lines from these 
transitions appear single but A 3980-9 shows two blended components. 
Therefore the 5D terms have all narrow structures. The width of A 4396-4 
enables an estimate to be made of the structure of 54°D, while from A 3980, 
the total width of 4d5D, can be estimated. The line A 3935 gives an 
approximate structure of the term 4d°D,. There is no check available. 


Fig. 1 
5165 5258 ; 5182 
F 
5p 
25 
40 
3 5 3P, 
F pas 3p A “18 
288: %e fe y Tr 
Vp Yo 
iS Ze 22 
5398, % 42 


em.—! x 10-3 
022 40 78° 103 


em,.—! x 10-8 


38 
| re # 
Cale. 60 © 38 0 0224 62647102 Deere | 


Vp 


in) [| 


48 0 violet 0 95 violet 0 38 93 violet 


§5. Discussion. 


Table ILI. shows the term hyperfine structures derived. A minus sign 
before the first interval implies that the whole structure is inverted, the 
highest F value lying deepest. 

The term interval factors derived from these structures are given in the 
last column of the same table. 

The 55°S,, 68°S, and 55S, terms are spherically symmetrical and are 
not therefore affected by any nuclear electric quadrupole moment which 
leads to deviation from the interval rule in the are spectrum of bromine 
for some lines. For these terms, the interval rule should be strictly 
obeyed unless perturbations are caused by configuration interaction. 
The structures of these terms as determined here show that the interval 
rule is obeyed within the limits of measurement, 
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Fig. 2. 
4316 4785 4704 
fe Ye 
5p 
2 3/2 c 
Yo 3h 
5p*P, 
Ie Vp 
1S3 
Sp eta apa itc aie 
5888, (09 
&6 “Ve : ee 
He e 
cm. x 10-3 | | 
Cale. 
0 153 = 260325 125 212 268 280 240 129 220 274 
em. x 10-3 
Obs. 
0 154 264323 213 270 0 130 219 274 
Violet violet violet 


TABLE IIT. 
Hyperfine Structures in the terms of BrI{—Intervals. 


Structure Interval factors 


Configuration Term Grete 108 cm 10 


4s? 43 5s Se 153 109 66 44. 

: 38, —38 pepe —14 

48? 43. 5p ePs 24 LSet. 5 

As 28 DA 112 8 

oP 2 ] 1 

2) Sof) ees Fae 6 

sp, 40s 16 

3P, 0 

4s? 43. 6s “eh 36 reg ag) 10 
4s? 43 4d 5D, Very small. 

5D, Very small. 

BD) wie 2 eee eae 9 26 

OE 10 6 4 

5D, 0 


The interval factor for the term 5s°S, is about four times that of the 
term 6s°S,. This is similar to what was found by Tolansky and Forester 
(1938) for the first spark spectrum of iodine. These terms arise from the 
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addition in parallel of an s-electron to the electron core of the basic ion 
term 4p°.4S5,. of the doubly ionized bromine atom, according to the 
coupling theory of White (1930). Thus from the smallness of the 
interval factor of the 6s°S, term, it follows that the coupling of the 4p* 
group is small. The term °S, is formed by the addition of an s-electron, 
anti-parallel to the 4p? group of BrIII. The interaction energy 1s 
proportional to cos JS which is equal to —1 in this case. Thus the 
interval factor for the 8, term should be negative and numerically smaller 
than the positive interval factor of the term °S,. The value of the 
interval factor for 8S, as found here (—14 cm.~!x 10-%) fits in the above 
theory. The magnitude, however, suggests that the coupling constant of 
the 4p electrons is about 4/7ths that of the 5s-electron. This is not as 
small as one would expect from the interval factor of the 6s5S, term. 
The interval factors of the 47°D terms are small and this suggests that the 
coupling of the 4p? group is not strong since a 4d electron is not expected 
to possess penetrating properties. In contrast to singly ionized iodine, 
the coupling of the 4p? electron group in BrII is weak and the interval 
factor of the 5s°S, term rather larger than would be expected under these 
circumstances. In BrI, the group 4p° has a large coupling constant. 
This may mean that there is a change in the nature of the electron coupling, 
in moving from the 4s?.4p* to 4s?.4p* configuration. This can also 
be inferred from the fact that contrary to expectations the singly ionized 
bromine atom shows structures which are no wider than those of the 
neutral bromine atom. 

No anomalies appear except for A 4356-9. It is a transition from the 
41°D, to the 5p°P, term. The structure of this line, which appears as 
a sharp single, does not fit in with the observed structures of the two 
terms involved. It is possible that the line is wrongly classified. 

Regarding the unclassified lines these can be arranged in groups’showing 
similar structures. The two lines A 4223-9 and A 3986-3 have similar 
structure suggesting that they involve the same common term. This is 
substantiated by the Zeeman splittings given by Lacroute. The lines 
A 4719 and A 4365 both show a triplet structure roughly of the same order 
while the three lines A 4297, 4179 and 3968 have similar patterns. From 
the data on Zeeman effect given by Lacroute, it seems that none of these 
lines can belong to the classified levels of the 4S system. 
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ABSTRACT, 


A brief review of published papers on the subject of electrical resistance 
at very low temperatures shows that several metals exhibit a feature 
which has not yet been explained theoretically, namely that the 
resistance goes through a minimum value as the temperature decreases, 
whereas present theories indicate that it should remain constant. An 
apparatus is described for measuring electrical resistance below 1° K. ; 
by means of its use, the effect was confirmed in gold and silver and also 
found in copper. There was no tendency for the curves to flatten off even 
at the lowest temperatures. The gold behaved in magneto-1esistance 
experiments as if it contained less than 1 part in 10° of ferromagnetic 
impurity. Slight anomalies in the shapes of some of the curves seem also 
to exist. 


$1. INTRODUCTION. 

Ar very low temperatures it is regarded as usual behaviour for the 
electrical resistance of a non-superconducting metal specimen to fall 
to a value which is constant, independent of temperature. The accepted 
explanation is that when the thermal vibrations of the lattice are 
negligible, the mean free path for electron scattering, which determines 
the electrical resistance, is constant with temperature, being in turn 
limited only by the imperfections in the crystal lattice. The imperfections 
include foreign impurity atoms and “ physical” impurities such as 
dislocations, a mosaic crystal structure or other departures from perfect 
regularity, which do not change with temperature. 

A considerable number of exceptions to this rule have however been 
noted. Among metals whose resistance begins to increase again as the 
temperature is decreased, the best known is gold which has been the 
subject of several investigations at liquid helium temperatures (Meissner 
and Voigt 1930, de Haas, de Boer and van den Berg 1934, de Haas and 
van den Berg 1936, 1937, Giauque, Stout and Clark 1937, de Haas, 
Casimir and van den Berg 1938, Stout and Barieau 1939, Nakhimovich 
1941, Garfunkel, Dunnington and Serin 1950). Of a number of specimens 
of silver which have been investigated, one was found to exhibit the same 
phenomenon (de Haas and van den Berg 1936). Similar behaviour was 
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noted in gallium by Keesom (1933), and in sodium, magnesium, aluminium, 
molybdenum, cobalt, and cerium by Meissner and Voigt (1930). Sodium 
had been previously found to exhibit an increase with decreasing 
temperature by Woltjer and Kamerlingh Onnes (1924); the specimen 
measured by Justi (1948) exhibited a decrease at low temperatures 
while MacDonald and Mendelssohn (1950) found ‘‘ normal” behaviour 
for a purer specimen. Keesom (1934) and more recently Garfunkel, 
Dunnington and Serin (1950) and MacDonald and Mendelssohn (1950) 
have confirmed the effect in magnesium. Keesom (1933) apparently 
observed a minimum for aluminium, but Boorse and Niewodniczanski 
(1936) and Garfunkel, Dunnington and Serin (1950) found no minima for 
their specimens of aluminium, while Steiner and Fiinfer (1936) found none 
for molybdenum. 

The phenomenon in gold is regarded as one not merely due to spurious 
causes ; the effects of varying the measuring current, of attaching the 
electrical leads to the specimens differently, and of varying the heat 
treatment of the metal were all eliminated by de Haas and van den Berg 
(1937). The phenomenon still occurred in specimens for which no trace of 
ferromagnetic impurity was quoted, but the possible effect of such 
impurities may nevertheless be important and is dealt with later in the 
present paper. 

Among suggestions as to the cause of the phenomenon there are those 
due to Justi (1942) that it is a size effect and Gorter (1938) that it may 
be a thermodynamic necessity—both of which have been criticized by 
MacDonald and Mendelssohn (1950)—while Lane (1949) has suggested 
that scattering at grain boundaries may play a part, though this 
mechanism would not seem to lead to a variation of resistance with 
temperature. The theory of Vonsovsky (1940) apparently refers to 
electric fields which are many orders of magnitude greater than those 
encountered experimentally. Sondheimer and Wilson (1947) have 
treated the problem of electrical resistance using a model more 
complicated than the free-electron model and have shown that at 
intermediate temperatures the resistance cannot be split up into two 
additive components (temperature-variable and temperature-independent). 
They assume from the outset, however, the existence of a temperature- 
independent residual resistance at very low temperatures. The 
interactions introduced by Frohlich (1950) into the free electron model 
to account for the phenomena of superconductivity seem to have no 
bearing on the present problem. 

It is an open question whether the rise in resistance continues as the 
temperature is decreased, particularly below 1° K. The only experiment 
so far performed in this region is that of Casimir, de Haas and van den 
Berg (1938) on gold. They showed that the rise did continue but the 
accuracy of their results was apparently impaired by poor thermal 
contact between the gold wire and the paramagnetic salt used to attain 
the low temperature. The present experiments are an attempt to repeat 
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this work—measurements being also taken throughout the liquid helium 
range, 4:2° K. to 1° K.—and to extend it to other metals ; this report 
presents preliminary results. 


§2. APPARATUS. 


A method has been described (Mendoza 1948) for establishing thermal 
contact between a metal specimen and a paramagnetic salt cooled by 
adiabatic demagnetization below 1° K. The metal specimen is joined to 
copper fins of comparatively large surface area embedded inside the pill 
of compressed powdered salt. When the salt is demagnetized, the metal 
cools. If power is then generated in the metal specimen—for instance, 
by the passage of a measuring current—a temperature difference is set 
up between the metal and the salt. Other things being fixed, this 
temperature difference increases rapidly with decreasing temperature ; 
its order of magnitude may be estimated, and it can be made small if 
the dissipating area of the fins is made large enough. At the same time, 
the provision inside the salt of layers of metal of comparatively high 
thermal conductivity is advantageous, reducing the inhomogeneities of 
temperature which might otherwise persist inside the salt and lead to 
errors of measurement. 

An apparatus embodying these data was constructed for resistance 
measurements ; the low-temperature part, which was immersed in a 
liquid helium bath, is shown in fig. 1. The salt A—iron or chrome alum— 
was cylindrical in external form, 12-5 mm. in diameter and about 30 mm. 
long, and contained four or more copper fins F, 5 mm. wide, 0-1 mm. thick, 
placed parallel to one another, running the whole length of the salt and 
projecting from the ends. The fins had a number of small holes cut in 
them to increase the mechanical strength of the block of salt; their 
total surface area was approximately 12 cm.”._ The assembly was made in 
one operation by pressing in a suitable die. M was a former, a hollow 
cylinder of copper, 10 mm. external diameter, 0-5 mm. thick and 16 mm. 
long. It had a spiral groove on its outside, deep enough to hold, without 
strain, the wire W whose resistance was to be measured, and the surface 
was insulated with a thin covering of a bakelite varnish. M was connected 
by a wire to one of the fins. With this construction consistent results 
were obtained. The ends of the specimen W were joined to short 
potential and current leads connected separately to the fins. The other 
ends of the fins were soldered to superconducting leads (of which one is 
shown at L), made most conveniently from 0:03 mm. diameter constantan 
wire, tin plated electrolytically then covered with a thin layer of soft 
solder. Such leads were found to be robust, flexible and easy to solder 
and to have a high thermal resistance. They passed from the fins, along 
the outside of the salt and joined on to platinum wires which passed 
through platinum.seals such as P into the helium bath, the length of lead 
between the salt and the platinum wire being about 5 cm. The magnetic 
susceptibility of the. salt was’measured by a static (Gouy) method, since 
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this ensured freedom from the eddy current heating in the fins and wire 
which might have occurred: had dynamic methods been used. A was 
suspended by a nylon thread from a quartz rod Q, attached in turn to 
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a Sucksmith balance (Sucksmith 1929, Jackson 1933) at the top of the 
cryostat. A number of flat spiral springs retained the quartz rod centrally 
inside the glass down-tube D; the lowest of these springs, 8, also served 


—,;, 


Resistance of Gold, Silver and Copper at Low Temperatures 295 


to maintain the top of the nylon thread at the temperature of the helium 
bath. The assembly was housed in a thin-walled calorimeter C, G being 
a copper-glass seal. Contact gas could be introduced into the calorimeter 
and pumped away again. 

The dimensions of the wire specimen W were chosen as follows. The 
smallest voltage that could be measured with the given potentiometer 
to the required degree of accuracy (50 nV. to 0-2 per cent in our case) 
was one datum, the other was the need to limit the power dissipation 
in the wire to a value determined by the area of the fins and the maximum 
tolerable uncertainty of temperature (in our case, the order of 10-1 erg/sec. 
giving a difference of less than 10-* deg. at 10-1 °K. and 3x 10-3 deg. at 
6x10 °K.; the heat conducted down the leads L also had to be 
dissipated but this was calculated to be less than 10-2 erg/sec. and was 
negligible). These conditions defined the resistance of W and the measuring 
current (0-25 ohm and 0-2 ma.), and, when using good conductors like 
gold, compelled the use of long thin wires. It was then necessary however 
to check that two further conditions were satisfied. The first was that 
the diameter so determined was sufficiently great, compared with the mean 
free path for electron scattering, for the effect on the electrical resistance 
of scattering by the surface of the wire to be negligible (in our case the 
diameter was usually at least 10 times the mean free path). The second 
condition concerned the temperature distribution along the length of 
the wire; had the wire been thermally isolated from its surroundings, 
except for ends maintained at the temperature of the fins, its temperature 
would have varied approximately parabolically along its length because 
of the generation of Joule heat at all points. In practice, an imperfect 
thermal contact probably existed with the former M, and the value of 
the temperature inhomogeneity calculated from estimates of the thermal 
conductivity of the wire represented an upper limit. Nevertheless the 
precaution was sometimes taken of connecting one or two points along BW 
to separate fins. (With sections 20 cm. long, the maximum inhomogeneity 
of temperature was calculated to be of the order of 3 x 10~° deg. at 1077 °K. 
and 5x10-% deg. at 610-2 °K. for most of our specimens.) A final 
‘restriction on the experimental conditions was to note that the warming 
time of the salt after demagnetization was sufficiently long compared 
with the thermal time-constant of the salt calculated from its known thermal 
diffusivity (Garrett 1950) and the linear dimension between fins. (This 
was easily satisfied, even for warming times as low as two hours.) 

The above conditions allowed measurements of known maximum 
uncertainty to be made below 1° K. For measurements on the same 
specimens at liquid helium temperatures, contact gas was introduced 
into the calorimeter and the same measuring power introduced negligible 
errors of temperature. Those specimens which were measured only in 
the helium range were immersed directly in liquid helium ; they were 
wound on a wax-coated mica cruciform, the wax being later dissolved 
away to leave the wire hanging freely. Again, during measurement the 
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rise of temperature of the specimen above that of the bath, due to the 
Kapitza sub-surface layer of poor thermal conductivity (Kapitza 1941) 
was kept to a low value with the small power dissipation. (In some 
measurements done in the past on short, fat specimens of high electrical 
conductivity, this condition may not always have been satisfied.) 


§3. ResuLts AND DIscUSSION FOR GOLD. 
Measurements in Zero Magnetic Field. 


The variations of resistance with temperature for three specimens are 
shown in fig. 2. All were from the same batch of material, the wires being 


Temperature (deg. absolute). 


drawn by Johnson Matthey from their Lab. No. 3459 ‘“ by methods 
where there was virtually no likelihood that they should pick up ferro- 
magnetic impurities’, the spectroscopic analysis of the wire when 
0-5 mm, diameter showing the presence of Fe and Ni (all sensitive lines) 
with Ag, Cu and Mg lines faintly visible and Sn and Na very faintly 
visible. The specimens were as follows: Curve I., diameter 0-025 mm. 
not annealed, measuring current 0:18 ma.; curve II. 0-025 mm. 
diameter, annealed in vacuo for two hours at 500° C., measuring current 


5-1 ma.; curve III., 0:10 mm. diameter, same heat treatment as IT., 
measuring current 20 ma, 
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It is seen that the resistance still continues to rise at the lowest 
temperatures reached. van der Leeden (1940) found a term proportional 
to the concentration of chemical impurity and varying as T~} to fit the 
Leiden results; our results do not fit any simple power law but vary 
more strongly with temperature than T-!. Our curves tend to rise more 
steeply with increasing minimum resistance, in agreement with the 
observations of previous workers. . 


Magnetoresistance Experiments. 


Magnetoresistive effects are conveniently shown on Kohler diagrams, 
where the fractional increase in resistance 4R/R is plotted as a function 
of H.Roog/R, where H is the applied field; curves for different 
temperatures are then usually superposable (see for example Blom (1950), 
where however the curve for gold is wrongly drawn). However, Stout 
and Barieau (1939) found that for a gold wire of high minimum resistance 
(R/Roo¢, ~ 650 x 10-4), which exhibited a marked rise of resistance in 
the helium range, the magnetoresistance curve for transverse fields at 
1-5° K. lay well below that at 4:2° K. The indication was that the effect 
would have been zero at approximately 1° K.; below this temperature 
it would presumably have changed sign and a negative magnetoresistance 
effect was in fact found at 1-5° K. in another specimen by Giauque and 
Stout (1938). Previously Meissner and Scheffers (1929) noted a negative 
effect at one temperature, 4:2° K., with a gold wire in small fields. The 
anomaly was attributed to ferromagnetic impurity in the wire and this 
was confirmed by Nakhimovich (1941) who made a careful study of the 
effect of adding iron impurity to the gold. He found that the deviations 
from Kohler’s rule, including a negative magnetoresistance effect, could 
be observed with gold wires in which of the order of 10° per cent of iron 
was dissolved ; the presence of the iron made the metal paramagnetic 
at low temperatures, and of high minimum resistance. The magneto- 
resistive change could be considered as consisting of two approximately 
additive components: a “normal” positive component, and an 
‘anomalous’ negative component, which was due to the presence of 
the iron ; the positive component could be increased with respect to the 
other by annealing. 

We have measured the transverse magnetoresistance effect in fields up 
to 10 kilogauss with gold wires from the same batch as those used for 
measurements in zero field. No anomalous effects were noted, curves 
from 14° K. down to 1-2° K. being shown in fig. 3; the facts that these 
are closely superposable and that plotted on a reduced Kohler diagram 
(Blom 1950) they join on to high-temperature curves for gold and silver 
seem, in view of the results of Nakhimovich, to indicate that the specimens 
were free from ferromagnetic impurity—even more free than N akhimovich’s 
purest specimens. It is possible to make a rough estimate of an upper 
limit of the concentration of iron in our specimens, by assuming first that 
Kohler’s rule holds for iron-free specimens, and secondly that the 
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percentage difference between the magnetoresistance effects at approxi- 
mately 4:2° K. and 1-:3° K. for specimens containing iron is proportional 
to the concentration. On the basis of the trend of the points in fig. 3, 
compared with the data for Nakhimovich’s specimen 13 and Stout and 
Barieaw’s specimen, 10-4 per cent is set as this upper limit. Though 
this does not demonstrate that the rise of resistances in zero field takes 
place in the absence of all ferromagnetic impurity, it would nevertheless 
seem surprising if 1 part in 10° of dissolved iron impurity were the cause 
of increases plotted in fig. 2. 


XiOer 


§4. REesuLTs AND DIscusSION FOR SILVER. 
Measurements in Zero Magnetic Field. 


The silver specimens were all 0:05 mm. diameter, drawn by Johnson 
Matthey from their Lab. No. 3406. Spectrographic analysis of the 
original silver showed the presence of Fe, Mg, Mn, Pb, Si, Cd (faintly 
visible) and Sn, Na, Cu (very faintly visible). Some resistance-temperature 
curves are shown in fig. 4. Details of the specimens are as follows : 
Curves I., II. and IIT., specimens not annealed ; measuring currents 14:5, 
2:5, and 0-18 ma. respectively. Another specimen Ag3 gave a curve which 
coincided within very small limits with curve III. between 4-2 and 1-23° K., 
(12 ma. measuring current.) It was then annealed for two hours at 450° (. 
in a vacuum of better than 10-° mm.; curve IV. was then obtained with 
14-5 ma. measuring current. It was then given three further heatings 
to 600° C. in helium gas; curve VI. was later observed using 14:7 ma. 
measuring current, This specimen was also used for magnetoresistance 
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experiments. Curve V. was for a specimen given the same heat treatment 
as that for curve IV., with 21 ma. measuring current. Curve VII. is 
taken for comparison from de Haas and van den Berg (1936) and is for 
a specimen of Hilger silver drawn down to 0-06 mm. diameter, annealed 
in vacuo at 450-500° C. for three hours ; measuring current not stated. 
The increase of resistance with decreasing temperature occurs with 
all specimens, implying minima in the curves above the helium range. 
The curves do not fit any simple power law. There is a general tendency 
for them to rise more sharply as the minimum value of the resistance 
increases, but this is not observed exactly. As remarked by de Haas 
and van den Berg (1937) concerning their work on gold, variations in 
the shape of the curves for unannealed specimens are not unexpected but 
there exists also a lack of homogeneity between the curves for the annealed 
specimens, which one might have expected all to have belonged to the 
same family. The anomalies observed by MacDonald and Mendelssohn 
(1950) in the resistance-temperature curves of potassium and cesium may 
be in some way analogous to these variations. It is also possible that 
the effects of dissolved gases, undetected spectroscopically, but capable 
of varying in amount and distribution on annealing, may be important. 


Magnetoresistance Expervments. 


The specimen measured by Stout and Barieau (1939) (which in zero 
field gave R/Rgo ¢ =194:2 x10, constant throughout the helium range) 
followed Kohler’s rule in its transverse magnetoresistive behaviour, their 
curves for. temperatures from 20-:3° K. down to 1-54° K. being closely 
superposable ; the dotted line in fig. 5 is for this specimen. 

Our specimen Ag3 however gave different results, similar in some but 
not all respects to the results for gold contaminated with iron. In fig. 5 
are Shown the curves in transverse fields at 4:2° K. and 1-23° K. for the 
specimen before annealing. After the first annealing however the effect 
plotted in this way was quite unchanged, though the resistance at any 
temperature had been greatly reduced. As mentioned above Nakhimovich 
noted a different result for gold, namely a tendency to reduce the departures 
from Kohler’s rule after annealing. It is conceivable that this difference 
of behaviour is in some way connected with the fact that iron, though 
it can dissolve in gold, is insoluble in silver (Constant 1945) so that the 
effects of annealing on the impurity are likely to be very different. 


§5. Resutts AnD Discuss1oN—COopPPER. 


Three samples of enamel-covered commercial copper wire were 
investigated, and the results are shown in fig. 6. All samples showed 
traces of Ca, Cd, Mg, Si, Pb, Mn and Zn on spectroscopic analysis, the 
total impurity content being estimated at less than 10-2 per cent. The 
wire for Curve I. was 0-035 mm. diameter, and the measuring current 
0-2 ma.; for Curve II. 0-05 mm., 0-18 ma.; for Curve III. 0-10 mm. 
4-5 ma, None was annealed, The comparative lowness of the semidnel 
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resistances was unexpected when Curves I. and IT. were measured, with 
the result that the measuring accuracy for both these curves was not 
high. The rise of resistance at the lowest temperatures is nevertheless 
outside experimental error. 

One further sample from a reel of copper wire, 0-035 mm. in diameter, 
of curious composition was also investigated. It contained 0-94 per cent 
of Sn, and also showed traces of Mn. It had a very high resistance at 
helium temperatures (Allen 1933), and the following values of R/Rgo ¢. 
were noted: 4:19°K., 0-4322; 2-50°K., 0-4324; 1-46°K., 0-4325; 
0-418° K. 0:4326; 0-067° K., 0-4327. 

Previous investigations of copper have given the following results. 
The readings of Meissner and Voigt (1930) show variations at liquid 
helium temperatures which are probably due to experimental error. 
de Haas and van den Berg (1934) found that the curves for several samples 
of pure copper were not parallel at liquid hydrogen temperatures ; 
van den Berg (1948) found that unannealed technical copper wires gave 
constant residual resistances at liquid helium temperatures. 

No magnetoresistance measurements have so far been taken on copper. 


§ 6. CONCLUSION. 


It appears from previously published work and from our own experiments 
that the rise of resistance at very low temperatures is a phenomenon 
occurring in a fairly large number of metals. It has only been noted so 
far in specimens which were far from being ideally pure, but in any case 
theory gives no indication as to the cause of the rise in connection either 
with pure or impure metals. It remains undecided; however, whether 
the rise is to be associated with the presence of small amounts (less than 
| part in 10° in our gold specimens) of ferromagnetic impurity. It also 
appears to be a fairly common occurrence that resistance-temperature 
curves for different specimens of a given metal do not always vary from 
one to another in a simple manner as the impurity is varied (see for 
example the curves for copper given by de Haas and van den Berg (1934) 
and our own curves for silver (above)). Presumably this means that 
different Nee may produce curves of the same minimum value of 
resistance but of different shapes. 
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XXXIV. On Pair Production by Fast Electrons. 
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SUMMARY. 

A study has been made of electron “ pairs ” directly produced by fast 
electrons in photographic emulsions. Evidence is presented to show that 
this process is of importance in the development of electromagnetic 
cascades present in the “ soft > component of the cosmic radiation. It is 
found that the variation with energy of the cross-section for the inter- 
action agrees with that predicted by theory. 


§ 1. INTRODUCTION. 


OBSERVATIONS on the creation of electron “ pairs” by fast electrons in 
photographic emulsions have been reported by Occhialini (1949) and by 
Bradt, Kaplon and Peters (1950). In this laboratory we have now 
searched, with immersion objectives, 2-8cm.* of electron-sensitive 
emulsion exposed to the cosmic radiation at an altitude of 22 km.; and 
we have found 80 examples of these so-called “trident ” or “ giraffe ~ 
events. Of these, the path-length in the emulsion, either of the “ primary ”’ | 
or of the three “ secondaries ”’, is in 50 cases of sufficient length to permit — 
a determination of the energy of the primary electron by the scattering — 
method. 

In favourable examples it is possible to make measurements of the 
energy of all the particles associated with a given “‘ trident’. In almost 
all such cases, we have found that the energy of the primary particle is 
equal to the sum of the energies of the out-going particles, within the limits 
of experimental error. Such a conservation of energy is confirmed by the © 
work of Dilworth et al. (1950), who have, in addition, found strong evidence — 
that charge is also conserved in the interaction. In those events in which 
the energy of the primary particle cannot be directly measured we can, 
in favourable cases, infer it from the total energy of the three secondaries. 
In determining the energies of particles we have adopted scattering 
constants and procedures recently developed by our colleagues in this 
laboratory and shortly to be published. 

The formation of secondary electrons as a result of interactions between 
fast charged particles and the Coulomb field of nuclei was suggested by 
Oppenheimer. Studies of theoretical aspects of such processes have been 
made by Furry and Carlson (1933) and many other authors. Bhabha 
(1935) has extended the treatment to include the effects of screening and 
gives an extensive bibliography to the earlier work. 


‘ 
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§ 2. IDENTITY OF THE TRIDENT PARTICLES. 


Two examples of tridents are shown in Pl. XIII. In A, two of the 
secondary particles originating at P are of relatively low energy and 
cannot have a rest-mass greater than 100m,. We may, therefore, with 
confidence, assume them to be electrons. Event B is typical of those in 
which the energy of the primary is shared almost equally among the three 
secondary particles. In such tridents, the origin is indicated by an 
apparent sudden increase in the grain density of the primary track, clearly 
seen at P. The collimation of the three tracks is so close that they are 
only just resolvable under high magnification at a point 300 microns from 
P; and they are not separated to the extent shown by event A until a 
point several millimetres from their origin has been reached. 

A number of observations give strong support for the view that at least 
most of the primary charged particles which produce pairs are electrons. 
It is well known that there are difficulties, in principle, in distinguishing 
between the tracks of particles, of charge |e| but of different rest-mass, 
when they are moving with great kinetic energy. Recent evidence from 
this laboratory strongly suggests, however, that the fast ‘ shower ” 
particles emerging from nuclear interactions are exclusively or almost 
exclusively, 7-mesons or protons (Camerini et al. 1950). We have never 
observed pair-production by such a particle. 

Further, the tracks of electrons of great energy can commonly be 
identified, if the length in the emulsion is greater than 2 cm., as a result 
of the reduction in energy of the particles through the creation of 
bremsstrahlung ; and we have frequently observed such effects in dealing 
with the primary particles of tridents. 

Finally, tridents are occasionally found to be produced by one of the 
components of an electron pair created by a photon ; and in favourable 
cases, when there is a considerable length of path in the emulsion, they 
may be observed to be associated with electron cascade phenomena. An 
example of such an event is shown in the line drawing reproduced in fig. ms 
In this case the primary particle of the trident is a fast recoil particle, 
almost certainly an electron, produced by an energetic heavy particle of 
charge |e|. The lateral scale in this diagram has been extended by a 
factor of 50, relative to the longitudinal, so that the characteristics of the 
different tracks may be displayed. The evidence for the electronic 
nature of the secondary particles is easier to obtain. They are frequently 
of such low energy that their rest-mass can be shown, as a result of scatter- 
ing measurements, to be less than 100m,. We therefore regard the evi- 
dence as very strong that both the primary and secondary particles of the 


trident events are electrons. 


§ 3. Tripents as A Disttvor PHysIcaL PROCESS. 
It has been suggested to us by Professor Peierls that the trident events 
might be due, not to the direct production of a pair of electrons by the 
charged primary particle, but to bremsstrahlung y-rays, which by chance, 
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have converted to electron pairs at points very close to the track of the 
parent electron (Hooper and King 1950). In order to examine this point 
we have made the following observations :—In the searched area of the 
plates we have found 150 pairs which evidently occur as a result of the 
cascade process in the emulsion-glass assembly. Events of this type may 
be seen in fig. 1. The perpendicular distance R, from the point of origin 
of the pair to the nearest parallel track, has been measured in each of these 
examples, and the distribution in R is shown in fig. 2. 


Fig. 1. 


Scems 


A reproduction of the tracks of a cascade shower plotted on a scale greatly 
extended laterally. AA represents the track of the heavy particle, and a 
trident is formed at P on the recoil electron track. The dashed lines indicate 
trident secondaries and the dot-dash lines the tracks of electrons from pairs. 


The magnitude and shape of this distribution at small values of R 
appears to be inconsistent with the view that the 80 tridents (R <0-25 p), 
are due to the same process. In such a case the distribution in the values 
of R would include those for the tridents (R <0-25 1) and would show a 
very sharp peak for very low values of R, followed by a slowly decreasing 
frequency from R=0-5 to 204. We have not calculated the distri- 
bution in the values of R to be expected from cascade theory, for it involves 
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extensive integrations, but it appears certain that it could not yield a 
distribution such as is obtained if the tridents are attributed to the 
conversion of photons. We therefore regard the observations as providing 
very strong evidence that most of the tridents are due to the direct pro- 


duction of pairs of electrons in the interactions of fast electrons with 
atomic nuclei. 


§ 4. THE Cross-SECTION FOR TRIDENT PRODUCTION. 

It appears to be possible, in principle, to determine the cross-section for 
the interaction leading to production of tridents, and its variation with 
energy, by measuring the total length in the searched area of the plates of 
electrons in different energy intervals ; and by comparing the results thus 
_ obtained with the corresponding numbers of observed tridents. Such an 
experiment demands, however, a very extensive investigation, for it 
depends on an identification of the electrons by the methods to which we 
have already referred. For this purpose only particles of great range in 
the emulsion can be employed and they are rarely observed. 


Fig. 2. 


PAIRS 


NUMBER of 


R in MICRONS 


Distribution in the values of the perpendicular distance, R, from the pair 
origin to the nearest parallel track. There are 14 events with R > 35 microns. 
The 80 tridents (R < 0-25 »), found in the same volume, would contribute a 
peak in the first interval. 


In these circumstances, we have determined the relative values of the 
cross-section as a function of the energy of the primary electrons. ‘This 
can be done since the distribution in energy of the electrons in the plates 
- employed has already been determined (Carlson et al. 1950). We assume 
the track length L—in the total volume of emulsion examined—due to 
electrons with energy between E and E-+dE to be represented by the 
equation 

We UN ADE OA OMA DE ie cel (1) 
where a is a constant, and y=—1-7. 

The distribution in the values of the energy of the primary particles of 
the 50 tridents on which it has been possible to make measurements, are 
shown in fig. 3. Fig. 4 shows the corresponding relative values of the 
cross-section for six energy intervals between 60 and 16,000 MeV. 
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It is possible to determine a mean, absolute, value of o over a very large 
energy interval, by considering those tridents which occur on one of the 
tracks of an electron pair created in the emulsion. These are then com- 
pared with the total track length of electrons, in the same energy interval. 
which originate in pair events in the emulsion. 
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The distribution, per unit energy interval, of the frequency of occurrence of 
tridents observed in Ilford G5 emulsions. 


Fig. 4. 


CROSS SECTION (ARBITRARY UNITS) 


O 10' ios 
Ey (Mev) } 


The observed variation with energy of the trident crogs-section.. The curve 
represents the theoretical variation normalized to the point at 1600 MeV. 


It may be shown that the expected number of tridents in the energy 
interval from E, to KE, is given by the relation 
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where a is defined in (1) and K is a constant depending on the absolute 
cross-section. Since we have observed only four tridents to originate on 
an arm of a pair, the value of 80 cm. obtained for the trident ‘‘ mean free 
path ” at E=10* MeV. is subject to a large statistical error, of the order of 
70 per cent. The observed value is, however, consistent with that cal- 
culated by Ravenhall; viz. 99cm. It will presumably be possible, for 
values of E up to 300 MeV., to determine an accurate absolute value for the 
trident cross-section in experiments with electrons from high energy 
accelerators. 
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Intrinsic Crystalline Structure and the Strength of Metals. 


By W. A. Woop, 


The Baillieu Laboratory, University of Melbourne*. 


[Received January 1, 1951.] 


In his paper ‘‘ On the Strength of Quasi-Isotropic Solids ’, Fiirth (1949) 
makes use of the idea that a solid contains an intrinsic block structure 
and derives a relation, similar to one obtained by Bragg, in which the 
strength varies inversely as the size of the blocks. Discussing the source 
of the structure, he suggests that the blocks might be accounted for by 
a recent theory due to Born, and goes on to imply that this theory would 
also explain the block structure in metals proposed by my colleagues and 
myself on the basis of X-ray diffraction work. This paper has been 
followed by a note by Paterson (1950) who queries these suggestions on 
the grounds that Born’s theory predicts a decrease in size of the blocks 
with increase of temperature. According to Paterson it ought then to 
follow from Firth’s hypotheses that a metal is stronger at higher 
temperatures, and that the X-ray diffraction lines from a deformed metal 
show greater broadening at higher temperatures of deformation. But, 
in practice, both these effects usually go the other way. 

It would be unfortunate, however, if the matter were allowed to rest 
at that, for the block theory of strength is the only one of any promise, 
being the only one that can be substantiated by direct observation. It 
might be worthwhile pointing out therefore that the arguments advanced 
‘by Paterson are less conclusive than they might seem at first sight. The 
essential point of the block theory should be that the strength, and the 
X-ray pattern, depend on the production in the structure of elements 
that are ‘dislocated’ or ‘‘disoriented”’; for only when they are 
dislocated or disoriented sufficiently can they provide sub-boundaries 
that are effective barriers to slip, or become incoherent X-ray reflectors. 
But there is no reason why these disoriented elements should always 
coincide with the ultimate blocks in any and every type of mechanical 
test ; they could be multiples of ultimate blocks ; and this aspect, as 
indicated below, is especially relevant in tests at higher temperatures. 

We think it is so important to distinguish between the actual 
‘ disoriented elements ” and the intrinsic ultimate blocks that it might 
be useful to illustrate the point by referring to the X-ray results. In 
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these, the fundamental observation is that the process of deformation 
reduces the grains of a metal to disoriented elements. If the deformation 
is drastic enough the elements become very small, but for a given metal 
they do not become smaller than a characteristic minimum size. Now it 
is evident that the breakdown could not be produced at all unless the 
grains in the first place contained enough lines of planes of weakness. 
Whether these faults form a static or dynamic pattern in the grain it 
is impossible to say. Nor, so far as the strength is concerned, does it 
matter ; for the only way of measuring the strength is by deforming the 
metal and as soon as the metal is deformed an observable disoriented 
substructure is produced at once. Therefore, for all practical purposes, 
we may speak of an intrinsic sub-structure, and this sub-structure has 
a lower limiting size. 

The next observation is that the size of the actual disoriented elements 
into which the grain breaks up is determined by the conditions of 
deformation, and in particular by the temperature and rate of strain. 
Thus, to quote from actual experiment: if aluminium of grain size 
10 cm. is strained at 300° C. at a rate of 0-1 per cent per hour, the 
resulting sub-structure is coarse; the size of the elements is about 
10°? cm., as against an ultimate crystallite or block size of less than 
10~-* cm. for deformation at room temperature. At first sight it might 
appear that at the higher temperature the pattern of faults in the grain 
before deformation is spaced more widely, so that the grains could break 
down to correspondingly coarse elements only. But this is not so ; for 
we can break the grains down to a finer sub-structure merely by increasing 
the rate of strain at the same temperature, and there is no doubt that by 
a sufficiently fast rate of strain we could break them down to the 
characteristic lower limiting size. 

Further, direct observations show clearly that the strength exhibited . 
by the metal under given conditions of deformation is primarily determined 
by the size of the actual disoriented elements produced by that 
deformation. Using the aluminium again as example, we find that during 
deformation at 300° C. at an extremely slow rate the grains do not break 
down at all; this is because the deformation presumably occurs mainly 
by viscous flow at the grain boundaries, but this aspect here is immaterial ; 
The relevant point is the parallel observation that in the absence of 
breakdown the metal exhibits little or no resistance to the applied strain ; 
its strength is negligible. This is one extreme. But if we increase the 
rate of strain at the same temperature we break down the grains and the 
sub-structure becomes increasingly finer as the rate of strain 1s greater ; 
and at the same time we find that the resistance or strength of the metal 
increases in proportion. Conversely, we can change the fine sub-structure 
back into a coarse sub-structure, still at the same temperature, by 
reducing the rate of strain again; then we find the strength pies 
by the metal is again reduced. Because of this eee agemt e 
sub-structure, temperature and strain-rate, 1f can be shown that a 
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specimen strained rapidly at an elevated temperature may be stronger 
than one strained slowly at a lower temperature. By experiments such 
as these it is not difficult to demonstrate that the strength follows the 
size of the sub-structure. 

These illustrations will therefore serve to indicate that the block theory 
of strength is on the right lines. It is true that in the form usually 
advanced by theoretical workers it is often open to criticism, of which the 
remarks by Paterson on .the Fiirth theory are an example. But the 
weakness of this theory as it stands (and of the corresponding Bragg 
theory) is not in the assumption of an intrinsic block structure, but in - 
failing to take into account the experimental observation that the 
operative elements are not always the ultimate blocks ; they may contain 
a group of blocks. It is also clear, however, that the block theory can 
easily be brought into line with experiment if so desired. 


Quadruple Moments and the Nuclear Shell Model. 


By B. TouscHEK, 
Department of Natural Philosophy, The University, Glasgow*. 


[Received January 22, 1951.] 


Ir has been pointed out by Hill (1949) that the quadrupole moments of 
atomic nuclei are in reasonable qualitative agreement with the nuclear 
shell model recently suggested by Haxel, Jensen and Siiss (1949, 1950) 
and M. G. Mayer (1949, 1950). We want to show that this model makes 
quite definite predictions about the quadrupole moments of a certain 
category of odd nuclei and that on the whole these predictions are in 
good agreement with experiment. . One should expect that the quadrupole 
moments of nuclei of the form (A—Z, Z)=(M,M-+1), where M is one 
of the magic numbers 8, 20, (28, 40), 50, 82, 126 should be determined by 
the odd proton alone, since the contribution to the quadrupole moment 
of the closed shells is zero. The periodic table shows two nuclei of this 
type : jgK3) and .,Biio%. The quadrupole moment from a proton in a state 
with spin j should be 


where the — sign corresponds to a nucleus (M, M+-1) as g2;,Bi the + sign 
to a nucleus (M,M—1) as 9 ,K. C=(1—(2A—Z)/A2) is a factor 
correcting for the motion of the mass centre and 7? is the average square 
distance of the odd (or missing) proton. It is obvious that the sign of the 
quadrupole moment predicted for K4! is wrong (as noticed by Hill), but 
the present experimental value is hardly more than an estimate. 


* Communicated by the Author. 
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Determining (r*)? for Bi one finds a value 7-4 10-! em. which is fairly 
near to the fast neutron radius (~9-10- cm.). <A further category of 
nuclei for which simple predictions can be made are odd proton nuclei 
with spin 3/2. Since the neutrons can only contribute very little to the 
quadrupole moment (via the centre of gravity motion, their contribution 
being Z/A? Xsmaller than that of a free proton) the quadrupole moment 
must stem from the protons. According to the shell theory the spin 
3/2 should only occur in connection with a shell of spin 3/2, so that an 
odd proton can only mean one or three protons in a 3/2 shell. The 
prediction of the shell theory is therefore 
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the + sign holding for the case of three odd protons. In the figure 


\/(r?) for such nuclei is plotted against A1’* and it is seen that the general 
trend of the curve is in good agreement with the curve (drawn out) 
‘representing the nuclear radii as measured by fast neutron experiments 
(Gamow and Critchfield 1949). It is worth noting that the positions of 
isotopes are according to their neutron numbers. Thus Cl?’ pip Te 
neutron number (20) lies lower than Cl®°, Ga‘ (40) lower than Ga® (38), 
and the positions of the Bromine isotopes may be due to the approach 
of the magic number 50. Bi with 126 neutrons lies comparatively 1G 
and the high quadrupole moment of As may be due to the small no. (2) 
of neutrons in the gg» shell. That the fluctuations in the case of 
quadrupole moments tend to be greater than those for the fast neutron 
radii seems to be a confirmation of a one particle picture, which will 
naturally lead to larger fluctuations than the statistical picture applicable 
to the determination of neutron radii, 
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That the odd neutron nuclei should have quadrupole moments 
comparable in size with the odd proton nuclei need not be a contradiction 
to the shell model, since the cancellation of the proton spins in pairs does not 
imply a simultaneous cancellation of the quadrupole moments. The 
check of these moments would require the determination of the con- 


tribution to the quadrupole moment of an arbitrary number of particles — 


in an unfinished shell. A few quadrupole moments (as those of Eu) are 
very large (about three times the size of the moments shown in the figure), 
and it is quite plausible to assume that they are the result of collaboration 
between three or more protons. This is further enhanced by the fact 
that none of the nuclei with ‘‘ abnormal” quadrupole moments has a 
spin <5/2. 
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The Velocity of Second Sound in Liquid Helium near the Absolute Zero. 


By J. C. Warp and J. WILKS, 
The Clarendon Laboratory, Oxford *. 
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RECENT measurements by Atkins and Osborne (1950) indicate that the 
velocity of second sound in liquid helium approaches the value 150 m./sec. 
at absolute zero, which is in agreement with Landau’s prediction of 
5/4/3 where C is the velocity of ordinary sound. Landau (1941) obtained 
this expression from the usual two fluid formula for second sound 
Cy=(TS*/C)(p,/p,,) where (p,/p,) was calculated in a rather obscure way 
from a model in which liquid helium at sufficiently low temperatures is 
represented by a phonon gas. However it is possible to derive the 
expression C//3 directly from the phonon gas model without having 
recourse to the two fluid theory at all. ) 
We follow Landau in representing liquid helium near the absolute zero 
by a phonon gas consisting of the elementary (Debye) excitations, for 
which E=pC where C is the velocity of ordinary sound and may be taken: 
as constant. We also assume that at these temperatures interactions in 
the helium are very small so that the collisions between phonons are 
elastic and conserve momentum. Actually Landau makes all these 
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assumptions in calculating (p,/p,). Suppose now that a compressional 
‘“‘sound”’ wave is propagated through this phonon gas, then there will be 
periodic variations of the phonon density. These periodic variations in 
the model will correspond to temperature variations in liquid helium, 
that is to second sound. 

Hence the velocity of second sound is the velocity of a sound wave in 
a phonon gas where the conservation laws of energy and momentum hold 
and also the relation E=pC. The expression for this velocity will be 
similar to that for the velocity of a sound wave in a photon gas for that 
too depends only on the two conservation laws and the relation K=pC;. 

The velocity of sound in a relativistic gas has been given by Hoffmann 
and Teller (1950) for the general case when a magnetic field is present. 
We may derive the same answer for the velocity in a phonon gas more 
directly by using a result due to Curtis (1950) that the velocity of sound 
in general relativity is \/(@P/0n) where 7 is the energy density. In a 
photon gas 7 is defined by YE=yC7 so by analogy in a phonon gas 


7nC=LH=v | Afexp (pe/kT)—1}"1. pe. dp... . (1) 
0 
Now the pressure of the gas is given by 


ae | (ve sin 6 cos 6 d0)(2p cos 6. Afexp . (pe/kT)—1 1-1 dp) 


=Joe | Ap{exp . (pe/kT)—1}* dp. Aaa 4: iinet?) 
0 


Comparing (1) and (2) 
P=4yC*, 
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This derivation gives no indication of the region in which we may 
expect this expression to be valid, Clearly it will fail if (a) the collisions 
between phonons do not satisfy the conservation laws of energy and 
momentum and (6) if the excitations of the liquid cannot be adequately 
represented by phonons. | All we can say is that if at a sufficiently low 
temperature the specific heat of the liquid varies as T® then we should 
expect the velocity of second sound to be C/4/3. It seems possible that 
this method could be used to predict the velocity at higher temperatures 
provided we knew a satisfactory expression for the other types” of 
excitation which may occur. Landau in his theory covers the region 
above 0-6° by introducing excitations called rotons but we can see no 
theoretical basis for his expression for them. ; 

As this model could also be used to represent a crystalline solid at low 
temperatures, it emphasizes a suggestion by Peshkov (1947) that second 
sound might be observed in crystals”. As he points out it will be necessary 


* It is interesting to note that Nernst (1917) suggested that in good thermal 


conductors at low temperatures heat might have sufficient inertia to give rise 
‘to an oscillatory discharge. 


ues: 
316 : Correspondence 


to experiment with a crystal in which the scattering of phonons by 
inhomogeneities and irregularities is a minimum. In a pure crystal 
however the essential criteria for second sound are probably that the 
number of elastic collisions between the phonons should be sufficient to 
ensure that the phonon gas is not in the Knudsen region while at the same 
time very few collisions of the umklapp type should occur. (This is 
important because in this kind of collision it is inadmissable. to apply 
the momentum conservation law (Peierls 1929)). Measurements of 
thermal conductivity by Berman et al. (1950) and Wilkinson and Wilks 
(1951) have shown that the number of umklapp type collisions rapidly 
decreases as the temperature falls and should be negligible in the region 
where the conductivity is limited by size effect. Therefore we intend to 
look for second sound in a corundum crystal at helium temperatures 
where umklapp collisions should be rare. 
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